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CARTA DE APRESENTAÇÃO 


Caro (a)s aluno (a)s, 

Sejam bem vindos a disciplina MATB97 - Equações Diferencias Ordinárias. Durante esse semes- 
tre, vocês terão oportunidade de estudar equações diferenciais e seus métodos de determinação 
de suas soluções. Esta discilina vai proporcionar uma experiência inicial com a teoria de equa- 


ções diferenciais ordinárias e modelagem matemática. 
E bem conhecido que muitos fenômenos que interessam às Engenharias e outras ciências po- 
dem ser estudadas através de modelos matemáticos nos quais aparecem de modo importante 


equações diferenciais ordinárias. 


Bom aprendizado e sucesso!! 
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UNIDADE 1 


EQUAÇÕES DIFERENCIAIS DE 
PRIMEIRA ORDEM 


Aula 1.1 


Introdução às Equações Diferenciais e 


algumas aplicações 


Por que as equações diferenciais são importantes? Onde elas aparecem ? Quando desejamos 
estudar Leis físicas que descrevem a natureza, usamos modelos matemáticos como aproxima- 
ção. Muitos destes modelos são relações que envolvem a taxa de variação de uma determinada 
grandeza física. Como foi visto no curso de Cálculo, a taxa de variação é uma derivada e as re- 
lações entre elas são equações. Sendo assim, estes fenômenos físicos são descritos por equações 
que envolvem derivadas que chamamos de equações diferenciais. 

A grosso modo, o que é uma equação diferencial? Sabemos que uma equação algébrica 
é uma equação que tem números como incógnitas. Numa equação diferencial, a variável 


incógnita é uma função y(x) e x é a sua variável independente. 


Definição 1.1.1. Uma equação que estabelece uma relação entre a variável independente x, a 
função incógnita y = f(x) e suas derivadas y',y”,...y!” se chama de equação diferencial. Pode 


ser escrita na forma 


PROC UG DU red vt) Ea): 


Exemplo 1.1.1. 
1. y = cosx, 


2. y+4x=0, 


MATB97 - Equações Diferenciais 


3. 22y"" +28y" = (2 +7)yº. 

A seguir, veremos alguns problemas que podem ser estudados usando as equações diferen- 
ciais. 
Modelo populacional 


Em um modelo populacional, dizemos que uma população cresce a uma taxa proporcional 
ao seu tamanho. Este modelo é razoável para população de bactérias, por exemplo. Que 


equação matemática representa este modelo? As variáveis do problema são: 


e téa variável independente e representa o tempo, 


e Péa incógnita e representa o número de indivíduos da população. 
Assim, o modelo populacional pode ser escrito da seguinte forma: 


e A taxa de crescimento populacional: AE 


e Se a taxa de crescimento da população é proporcional ao tamanho da população, então 


dP 


GE kP, 


onde k é a constante de proporcionalidade. 


Esta equação é uma equação diferencial que envolve a função incógnita P(t) e sua derivada. 


Note que P(t) = Ce“ é uma função que satisfaz 


dP 


Re kt = 
7 = MCe!) =kP. 


Sendo assim, esta função é uma solução da equação diferencial. Veremos estes conceitos com 
mais detalhes a seguir. 


O que significa a constante C? Note que, quando t = 0, 
P(0) =€. 


Podemos dizer então que C é a população inicial ou “condição inicial” do problema. 
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Lei de Hooke: um modelo para o movimento da mola 


Consideremos o movimento de um objeto com massa m na extremidade de uma mola 


vertical como mostrado na figura at 


[rara E 


2 ” ” 


| E 


A Lei de Hooke diz que se a mola for esticada (ou comprimida) x unidades a partir da 


Posição 


de equilíbrio 


sua posição de repouso, então a mola exerce uma força sobre o corpo que é proporcional ao 
deslocamento x: 


força elástica = —kx, 


onde k é a constante positiva da mola e x = x(t). 
Ignorando qualquer resistência do ar ou do atrito e sabendo que F = ma, podemos formular 


a equação da lei de Hooke: 


1.1.1 Classificação 


As equações diferenciais podem ser classificadas levando-se em conta as seguintes caracte- 


rísticas: 
1. O número de variáveis independentes da função incógnita; 
2. Onúmero de funções incógnitas; 
3. A estrutura da equação e 


4. A ordem da equação. 


MATB97 - Equações Diferenciais 








Exemplo 1.1.2. 





1 po + 20 = V(t), (modelo de circuito RC) é uma equação diferencial ordinária. Função 
incógnita Q = O(t). 

= a + a = 0. (Equação de Laplace) é uma equação diferencial parcial. Função incógnita 
u E Ra 


3. us + uu = O é uma equação diferencial parcial. Função incógnita u = u(x,t). 


Quanto ao número de incógnitas, teremos a seguinte classificação: 





Exemplo 1.1.3. Os exemplos dados em 1.1.2 são equações diferenciais com apenas uma função 
incógnita: 1- Q(t), 2— u(x, y) e3— u(x, t). O seguinte exemplo é um sistema que envolve as duas 
funções incógnitas x(t) e y(t) 

x (b) = y(b) 

vb) = (8). 





Joseph N. A. Yartey & Simone S. Ribeiro 


Exemplo 1.1.4. 


1. y=cosx => primeira ordem. 


2.y'+4x=0, => segunda ordem. 


3. 02y" +2ºy" =(02+79)y = terceira ordem. 
y y y 





1.1.2 Soluções de Equações Diferenciais Ordinárias 


Considere a EDO de ordem n 


EU sea VM) =0 (1) 





2 
Exemplo 1.1.5. Verifique que d(x) = z + -, sendo c qualquer constante real, é uma solução 


da equação diferencial xy + y =x? para todo o x E R. 
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Solução: A EDO pode ser escrita na forma 


Flx,yy)=xy +y-x2=0 


2 
2 
. Se P(x) = s + -, então a sua derivada é q'(x) = E — - Substituindo na equação, temos que: 
Fe 6660) = dE-S) (5 +L)-x 
MARA E 5 uRA Dee 


2 
é 
Isto mostra que a função (x) = CRE satisfaz a equação diferencial xy + y = xº, sendo assim 


uma solução da mesma. 
Uma solução também pode aparecer na forma implícita, isto é, como uma função da forma 


H(x,y) = 0. 


Definição 1.1.7. A relação H(x,y) = O é chamado uma solução ímplicita da EDO (1) se esta 
relação dá origem a pelo menos uma função de valores reais p(x) definida no intervalo I, tal 


que p(x) é uma solução explicita de (1) em. 


Veja o exemplo a seguir. 


Exemplo 1.1.6. Mostre que a relação x? + y? — c = 0, sendo > 1 constante, é uma solução na 


forma implícita da equação diferencial F(x,y,y) = yy' + x = 0 no intervalo (—1,1). 


Solução: A relação x? + y? — c = 0 produza as funções 
h(o)=y=NVc-x e kx)=y=-Vc-x? 


ambas definidas em (—1,1). 


—x 
Como h'(x) = ——, então 
c—x2 





F(x, h(x), h'(x)) = (Vc —2 ( 


) +x=0 eé definida em (-1,1). 
c—x? 
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Isto mostra que a função x? +y? -c=0, c > 1, satisfaz a equação diferencial yy' = —x, sendo 


assim uma solução da mesma. 


Exemplo 1.1.7. A relação x? + yº + 1 = 0 NÃO é uma solução na forma implícita da equação 


diferencial F(x,y,y') = yy +x =0 no intervalo (-1,1). 


Solução: A relação não produza uma função de valores reais no intervalo, resolvendo por y 
temos p(x) = y = V-4— x2 e este função não é definida pois —4 — x? < 0. 
Entretanto quando derivamos a relação implicitamente, temos que 2x +2yy' =00oux+yy' =0, 


que é a equação diferencial que queremos resolver. 


Observação. Vimos nos Exemplos 1.1.5 e 1.1.6 que uma equação diferencial PODE ter infinitas 
soluções. Se uma solução envolve uma constante real arbitrária c, como descrita nos Exemplos 1.1.5 e 
1.1.6, então ela é chamada de solução geral de uma equação diferencial. Se escolhermos uma constante 


específica c, encontramos uma solução particular do problema. 


Exemplo 1.1.8. Mostre que y = 2cx? + c?, sendo c uma constante, é uma solução geral de 


dy 3dy 2 
(28) OX q ARE y 


e ache uma solução particular que satisfaça a condição y = —1 quando x = 1. 


d 
Solução: Se y = 20x? + c?, então = 4cx. Portanto 


dy Ê dy 
3 2 3 
] + 8x E (4cx)* + 8xº(4cx) 


16c?x? + 32cx* 


16x? (20x? + 2) = 16229. 
a — 
y 


Portanto y = 2cx2 + c? é uma solução da equação diferencial. 


Se y = —1 quando x = 1, temos 
-l=2+0 >cC+%+1=05(c+1)/=0>c=+1. 


Portanto a solução particular é y = 1 — 2x2. 
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1.1.3 Uma família a n— parâmetros de soluções 
A solução geral da EDO de ordem n 
POC AU ter uv) =0 


é uma função que possui n constantes c1,c>,:-- , cn, que chamamos de parâmetros. Portanto a 
solução geral da EDO de ordem n é uma família a n-parâmetros de soluções (na forma explicita 


ou implícita) que contém todas as soluções possíveis num intervalo T. 


1) 


X 


Exemplo 1.1.9. Do Exemplo 1.1.5 temos que y = 2 
soluções de xy' + y = x2. 


é uma família a 1-parâmetro de 


“10 


Exemplo 1.1.10. Mostre que y = cj" + c,e 2 é uma família a 2-parâmetros de soluções de 


d 
e ache a solução particular que satisfaça a condição y = — = 1 quando x = 0. 


d2 


2 u—— 
dx? 


d 
Solução: Se y=cie" + c2e **, então a =q"-20e*?” e = qe + 40,e”. Portanto 
E? 


dy dy 
de “a É 


(ce + 40909) + (cqe” — 20,909) — Acre” + co 9) 


e +40e 7 +ce Ie” Ie — 2030 ** 


=0 


Portanto y = cj" + c2e * é uma família a 2-parâmetros de soluções da equação diferencial. 


d 
Se y = — = 1 quando x = 0, temos 


I 

pa 
PE 

pa 
De” 


Cc +c2 


II 
fes) 
— 
NO 
a 


Cj — 209 


Resolvendo temos cy = 1 e co = 0. Portanto a solução particular é y = e”. 
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1.1.4 Determinando EDO associada a uma família a n-parâmetros 
de curvas 


Seja uma família de curvas a n— parâmetro. Podemos considerar a família como a solução 
geral de uma equação diferencial de ordem n. Vamos tentar encontrar a equação diferencial 
correspondente. A equação diferencial da família, não deve ter parâmetros (constantes arbi- 
trárias) e sua ordem deve ser igual ao número dos parâmetros da família. A estratégia geral 
para encontrar a equação diferencial de uma dada família a n— parâmetros é derivando n vezes 
a equação. Isto produziria um sistema de n + 1 equações que podem ser usadas para eliminar 


os parâmetros. 


Exemplo 1.1.11. Dado a curva y? = cx + 3 a 1-parâmetro, determine a equação diferencial da 


família. 


Solução: Temos que: 


a =00+3 (1). 


Uma vez que temos uma família a 1-parâmetro, estamos procurando uma equação de primeira 


ordem. Derivando a equação temos 


2yy = 3cx (2). 
Encontrando o valor de c em (2), temos: 
2yy” 
— E (3 
= (3) 


Substituindo (3) em (1), temos: 


" 


2yy 2xyy” 
2 = [= Ig = 
y = 332 ) +3 3 +3 


po 3-9 


2Xy 





que é a equação diferencial da família. 


Exemplo 1.1.12. Determine a equação diferencial da família de curvas a 2-parâmetros y = 


— + Co. 
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cq e 
) X 
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(1) 


Uma vez que temos uma família a 2-parâmetros, estamos procurando uma equação de segunda 


ordem. Derivando a equação 2 vezes e encontrando o valor de cy, temos: 








Die Cc = 2 , 
y = ES DD (=X y 
” 20 . Ru 
DU e e 
Igualando (2) e (3), temos: 
210 Rus 
us > 


Simplificando, obtemos y” + 2x2y' = 0, que é a equação diferencial da família. 


d” 
1.1.5 Solução da EDO do tipo — = f(x) 





Lema 1.1.1. Se uma EDO de ordem n é da forma 


d'y 
dx" = f(x), 





então sua solução geral é obtida de uma maneira direta por integrações (sucessivas). 


Exemplo 1.1.13. Ache a solução geral da equação diferencial 


dy +2 

dg — C08X+2x. 
Solução: Seja 

dy 

— = cosx +21. 


dx 


Integrando, temos: 


E 
Il 


[cos x+ 20 ax 


> y = senx+xi+c. 
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Exemplo 1.1.14. Ache a solução geral da equação diferencial 


d2y Bd 

de rc 
Solução: Seja 

d2y ge a 

de rc 


Integrando, temos: 


d 
= = [ese =9 9 de=5 tr) +a. 


Integrando novamente, temos: 
y= fot +ea)dx=% +Inx+cax+c. 


Portanto, a solução geral éy= 27º +Inx+ ax +. 
ção 8 y 


1.1.6 Problema de Valor Inicial da EDO de ordem n 
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Exemplo 1.1.15. Resolva o PVI 


y”" -e + cosx=0 


(0) =2 

y(0) =3 

yv (0) =83. 
Solução: Seja 

dE a 

o =€ — COS J. 


Integrando, temos: 


dy : x 
q2= |(€- cosx)dx=e- senx+c. 
dx 


Integrando novamente, temos: 


dy 


q (*—- senx+c)dx=e + cosx+cx+c. 
x 


Integrando novamente, temos: 


1 
= + cosx+ax+o)dx=e+ senx+=X +cox+cs. 
4 2 


1 
Portanto a solução geraléy = e" + sen x+ 5010 + Cox + ca. 


Sey=2,y =3ey” =3 quando x = 0, temos 


l+c3=2 > c=1 
1+1+0=3 > 0=1 


1l++0,=3 > qo=2. 
Portanto, a solução particular éy= e + senx+1 +x +12. 


1.1.7 Exercícios 


1. Equações diferenciais ordinárias e parciais: Se a função incógnita da equação diferencial 
depende de apenas uma variável, então temos um exemplo de equação diferencial ordi- 
nária. Caso contrário, teremos uma equação diferencial parcial. Verifique se a equação 


diferencial a seguir é ordinária ou parcial e se a função dada é uma solução da equação: 
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(a) Uxx + Uyy = 0; u(x,y) = In(x? + y?). 
(b) y” -y=0; y(t) = é. 

(c) y" —-y=0; y(t) = cosht. 

(d) uxx + Uyy = 0, u(x,y) = cos xcosh y. 
(O) ty -y=Piy=3t+Ê 

1/20/40 


(£) Cu =UGu= (7/t) 


t 
(9) y —-2ty=1;y= e” á eds + e”. 
0 


2. Determine os valores de r para os quais a equação diferencial dada tem solução da forma 


=e. 
(a) y +2y=0. Resp: r = —2. 
(b) y'+y —-6y=0. Resp: r=2,-3. 
(c) y' -y=0. Resp: r= «1. 


(d) VW" => 3y" A 2y' = 0. Resp: r=0,1,2 


3. Determine os valores de r para os quais a equação diferencial dada tem solução da forma 


yv=t,parat> 0. 


(a) 2y" +4ty +2y=0. Resp:r=-1,-2. 


(b) Py” —-4ty +4y=0. Resp:r=lour=4. 
y Y y Pp 


- Equações lineares e não lineares. A equação diferencial ordinária 


Eli as VM) =0 


é linear se F é uma função linear nas variáveis y,1/,..., y!?. Esta definição é equivalente à 
de uma transformação linear vista em Algebra Linear. Uma definição similar se aplica às 
equações diferenciais parciais. Desta maneira, a forma geral de uma equação diferencial 


linear de ordem n é: 


ay? + a(Dy DD +. + ay = (0. 


Uma equação que não respeita a condição acima é não linear. Para cada uma das seguintes 


equações diferenciais abaixo, diga a sua ordem se é linear ou não linear. Explique. 
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(a) 2y” + ty +2y= sen t. Resp: segunda ordem linear. 

b) y"+y”" +y'+y +y=1. Resp: quarta ordem linear. 

(c) y" + sen (t+ y) = sent. Resp: segunda ordem não linear. 
(d) y' + tp = 0. Resp: primeira ordem não linear. 


tr 


+ ty +( cos 2Dy =p. Resp: terceira ordem linear. 


(e) y 
5. Mostre que y = a cos (mx + b) é uma solução da equação diferencial y” + m2y = 0. 


A dê. Da 
6. Mostre que v = B+ E é uma solução da equação diferencial - T=0 


de ads 


2 


7. Mostre que Ax?+B1? = 1 é uma solução da equação diferencial x Ga + dy E — dy =0 
, d y s duas Ya? "ax Vga 


8. Determine a equação diferencial da seguintes famílias de curvas. 
(a) y = e”*, onde m é constante arbitrário. 
(b) y = Ae” + Be”, onde A e B são constantes não iguais. 
(c) y = k arcsenx, onde k é um parâmetro 
B 
(d) y = Ad + E 
(e) y=e(Acos x+B sen x) 
(f) y=A+Be + Ce 
(g) xy= A + Be +xº 
(h) y-2ay+2 =0? 
9. Determine a equação diferencial de todas as retas que passam pelo origem. 
10. Determine a equação diferencial de todas as circunferências de raio R. 


11. Determine a equação diferencial de todas as retas que dista 1 unidade da origem. 


12. Resolva PVI 


dy dy a 
EU cóna 2 62-1=0 
3 9x = 12 (b) EE x 
(a) y(o) = 1 yvl)=5 
y(0) =2 
y"(0) = 4 
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d2y dty 


E o SO cos x 
(c) y(O) =1 (d) y(O) =0 
v'(0) = 4 y'(0) = y”(0) = —1 


y"(0) — 7 
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Aula 1.2 


Equação Diferencial de Primeira Ordem 


Definição 1.2.1 (Equação Diferencial de Primeira Ordem). Uma equação diferencial de 
primeira ordem contém somente a primeira derivada y' = = possivelmente y e alguma 


função de x, isto é, é uma equação do tipo 


Piva p=) ou pn) (IZOIA 


Exemplo 1.2.1. 
Te. E SOS; 


2. xy +4xy=0. 


1.2.1 Interpretação geometrica da EDO de 1º ordem 
d 
A EDO de 1º ordem — = f(x,y) associa a cada ponto (xo, yo) do domínio D € R? uma 
d 
direção m = N = f(xo,y0). A direção em cada ponto (xo, yo) é a inclinação da reta tangente 
à curva, com equação f(x,y) = c, passando pelo ponto (x9,y). O domínio D com a direção 
em cada ponto é chamado de um campo de vetores. Resolver a equação diferencial significa 
determinar curvas cujas retas tangentes no ponto (xo, 0) têm inclinação m = — = f(xo, Vo). 
dy 


Exemplo 1.2.2. A EDO + = y-— x forma um campo de vetores em R?, como mostrado na 


dx 


figura (a) abaixo. Cada curva y = x +1 + ce”, onde c é um constante arbitrário, possui uma 
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reta tangente no ponto (x, y) com inclinação m = y— x. A Figura (b) abaixo mostra uma solução 


; 2 
particular nassando pelo ponto (n. =. 





1.2.2 Problema de Valor Inicial da EDO de 1º ordem 





Exemplo 1.2.3. Podemos reescrever o Exemplo 1.2.2 como um problema de valor inicial 


dy 
de ?* 


nO) = 5 
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1 
cuja solução éy = x +1 — 3º 
Queremos analisar as seguintes questões: 
1. Sob que condições podemos garantir que existe uma solução do PVI (1.2.2)? 


2. Sob que condições podemos garantir que a solução de ( 1.2.2) é única? 


Teorema 1.2.1. (Existência ) 


Suponha que f(x, y) é continua num retângulo R = (x,y): xo-h <x<xo+h, yo-k<y< yotk) 


centrado em (xo, yo). Então existe um número h (possivelmente menor que h ) tal que uma solução 


y = Q(x) de (1.2.2) é definida no intervalo (xo — hm, xo + m). 





Demonstração. A demonstração desta será omitida, mas pode ser vista em [1]. 


Teorema 1.2.2. (Unicidade ) 


o 
Suponha que f(x,y) e sua derivada parcial io 4) sejam continuas num retângulo R como no 


Ox 


Teorema 1.2.1. Então existe um número h> (possivelmente menor que hy ) tal que uma solução 


y = q(x) de (1.2.2), cuja existência é garantida pelo Teorema 1.2.1, é a única solução de (1.2.2) 


definida no intervalo (xo — ho, xo + ho). 





Demonstração. A demonstração desta será omitida, mas pode ser vista em [1]. 


Exemplo 1.2.4. Considere o PVI 


dy 
(y + Ds = cosx 
(0) =2 


Verifique que a existência e a unicidade de soluções estão garantidas. 


Solução: Colocando a equação diferencial dada na forma exigida pelos Teoremas 1.2.1e 1.2.2 


temos: 
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o = 
Note que a função f e a sua derivada ef a — 
oy  (y+1) 





são continuas exceto quando y = —1. Então 


9 
Poe E são continuas num retângulo centrado em (0,2) que não corta a reta y = —1. O Teorema 
1.2.2 garanta que existe uma solução y = (x) em algum intervalo centrado em x = 0 e que esta 
solução é única em algum intervalo (possivelmente menor) centrado em x = 0. 


Provaremos mais adiante que a solução é p(x) = —-1 + V9+2. sen 1. 


Exemplo 1.2.5. Considere o PVI 


dy x 
dx y 
Y(1)=0 
Verifique que a existência e a unicidade de soluções estão garantidas. 
E E E O o o a se oREI ane E 
Solução: As funções f(x,y) = —— e o = E não são continuas em (1,0) (pois não são definidas 


em (1,0)). Então os Teoremas 1.2.1 e 1.2.2 podem ser usados para estudar a existência e/ou 


unicidade de uma solução para o PVI. 


dy X 
Observação. A relação x? + y = 1 satisfaz o PVI dx y , mas não é uma solução no sentido 
v) =0 


do Teorema 1.2.1. Recorde que uma solução deve ser definida num intervalo que contém o coordenada 
xo do valor inicial mas , y é definida implicitamente em x = 1 mas não num intervalo em torno de x = 1. 


De fato o PVI possui 2 soluções qi(x) = VI — x2 e qa(x) =— VI — 2x2, 


1.2.3 Exercícios 
1. Seja a equação diferencial x” = (x— 1) cos t e suponha que x é uma solução tal que x(1) = 1. 


Mostre que x(t) = 1 para todo t. 


2. Determine uma região no plano xy na qual a equação diferencial dada tenha uma única 


solução cujo gráfico passe pelo ponto (xo, yo) nessa região. 


d Eu = 2 
(a) us (co) (3y—- y) dy=(1+2xº) dx 
dy 
(6) qe = XY (d) (y- 9) dy = (y +) dx 
3. Considere o PVI : 
Ab io 
dx *y 
y(O) = O 
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1 
Mostre que y = 0e y = TE” são soluções. Será que isto contradiz o Teorema 1.2.2? 


4. Verifique a unicidade ou não da solução do seguintes PVI 


dy 444 dy 1/2 
o (E! º o [a 
v(O) = 7 Y(0)=0 
CU maio AR 17 
(b) | de A = (d) | E 
y2) =1 v(O) = 0 
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Aula 1.3 


Equações Separáveis 





Método de solução 


Para resolver a equação ( 1.3.1), integramos ambos os lados com relação a x: 


d 
[ewgar= | fooarso 


d 
Do Cálculo, = dx = dy. Substituindo na expressão acima, chegamos ao método de solução: 


[eanay= [ rosário (132) 
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Supondo que f e g são funções contínuas, as integrais existem. Calculando as integrais de 


ambos os lados, chegamos à solução geral da equação. 
Exemplo 1.3.1. Ache a solução geral da equação diferencial 9y y' + 4x = 0. 
Solução: 

e Primeiramente, separamos as variáveis: 9y dy = —4x dx. 


e Integramos ambos os lados da equação acima em relação a x: 


9% 4 
fovar- [axar > — = te > 42º +9y? =, q=c, 
isro é, a solução geral da equação diferencial é 4x? + 9y? = c1. 


Exemplo 1.3.2. Resolva a equação y =1+1y. 


Solução: 





e Primeiramente, separamos as variáveis: T = dx. 


+y2 


e Integramos ambos os lados da equação acima em relação a x: 





d 
[Jo > arctgy=x+c > y=tg(x+c) 


Assim, y = tg (x +) é a solução geral da equação diferencial. 


Exemplo 1.3.3 (Decaimento ou crescimento exponencial). O decaimento ou crescimento 
exponencial é modelado pela equação diferencial y = ky. Se a constante real k > 0, teremos 


crescimento e se k < 0, teremos decaimento. Resolva esta equação separável. 


Solução: 
e E 
e Primeiramente, separamos as variáveis: — = kdx. 


e Integramos ambos os lados da equação acima em relação a x: 
dy (kx+c1) ci pkx kex 
E kdx > Inly=kx+og > íy=0"U=edle*=ce 


E y = e“tc1) é a solução geral. 
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Exemplo 1.3.4. Resolva o problema de valor inicial 


Solução: 


dy dx 
e Primeiramente, separamos as variáveis: — = -—, 
y x 


e Integramos ambos os lados da equação acima para encontrar a solução geral: 


d 
[5-[5 > Inly=-Inbd+o => Epis dE > v=c/x, 


brl 


onde c = eº. 
e Substituímos a condição inicial (quando x = 1, y = 1) na solução geral: 
l=ci > c=1 > vy=1/x 


é a solução do PVI dado. 


1.3.1 Equações da forma y' = G(ax + by +) 


Se a equação diferencial é dada na forma 


a + by + 
e (ax + by + c), 


então podemos usar uma mudança de variáveis u = ax+by+c para transformá-la numa equação 


de variáveis separáveis. Observe que 


=ax+by+c => tie PA E i(Se a) 
u=ax yY+C de EE ou ARS a). 
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que se reduz a uma equação diferencial separável: 


du 


RG 


Exemplo 1.3.5. Resolva 


d 
a =(2x+3y-1)+4(2x+3y- 1) 2-2. 





dx 
Solução: 
-. du dy Red ã 
Seja u = 2x +3y-— 1 então o 2 + q Substituindo na equação temos: 
du a ut+4 
Ee u+4u “= a 
Separando as variáveis, temos: 
3 
- du = dx 
u 


Portanto 





E du= | dx > ufa aa 
us +40 Tia ia 


S mrpd =p 


Logo a solução geral é (2x +3y-— 1) +4 = eftte 


1.3.2 Exercícios 


1. Use o método de separação de variáveis para resolver cada uma das equações diferenciais 





abaixo: 
po Ly 
tango 1+x' 
poe 
EA de 


(O) y = ++. 
(d) 4 -=2xy=%. 
(e) y =x2/y. Resp: 3yº - 2% =c. 


(f) y +y senx=0. Resp: y! + cosx=csey+0etambém y =0 para todo o x. 
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amo eg o UA e y 
(g) e Er x +27 -eb=cy+re £o. 
dy 
(h) di 2x dE ony= 1 
dy 
(1) ae Y =y 
Que 
(6) dy — 10lnx 
dx xy-xy? 
dy x+2y 
(1) ax *€ 


2. Ache a solução dos seguintes problemas de valor inicial. 


1Ey 





(a) y = YO) =1. 
des 
() y = Y0)=1 


(O) y = +91 +23), y(0) = 1. 

(d) y' — 2xy =x, y(0) = 

(e) y =y(1-22), y(0) = —1/6. Resp: y=1/(2 —x— 6). 
(O) y=(1-223)/y u(1)=-2. Resp: y=—V2x-202 +44. 
(g) xdx + ye “dy = 0, y(0) = 1. Resp: y = [2(1 - e — 1]!2 
(h) dr/dO = aa r(1) =2. Resp: r=2/(1-2m 0). 


v+2y , Y(O) = —2. Resp: y = —-2In(1 + 2) + aJ'2, 


)y = 

3º tg ue +(1-e) sec2ydy =0 
y(In2) = 
E: a 


(1) = FR 
v(O) = 2 


3. Resolva o problema de valor inicial y' = 2y? + xy?, y(0) = 1 e determine onde a solução 
atinge seu valor mínimo. Sugestão: Use o teste da primeira e da segunda derivadas para 


determinar o ponto de mínimo. Resp: y = —1/((x2/2)+2x— 1);x = —2. 
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4. Resolva o problema de valor inicial y' = 2 cos 2x/(3 + 2y), y(0) = —1 e determine onde a 


3 
solução atinge seu valor máximo. Resp: y = Et vsen2x+1/4,x=7n/4. 


5. Resolva a equação 
dy ay+b 


dx  cy+d' 


onde a, b, ce d são constantes. 


Resp: y="y+5 mkay+b+ka 40,09 +40. 





6. Resolva as seguintes equações diferenciais usando o método de separação de variáveis, 


após fazer uma mudança de variáveis apropriada: 


(a) y =(x+y)?. Resp: y=-x+ tg(x+0). 
(b) y =(8x+2y+1)2. Resp: 8x +2y+1=2tg (4x+0). 
(c) (2x + 3y — 1)dx + (4x + 6y — 5)dy = 0. Resp: x+2y+3In2x+3y —- 7] =. 


(d) (2x — y)dx + (4x — 2y + 3)dy = 0. Resp: 5x + 10y + c = 3In[10x — 5y + 6]. 
Qy-3x+73+3 


2 
e(u+x+1? 


(e) y = 


(D y RES ES 
(g) y=2-342x-y+7 


7. Suponha que o valor y foi investido numa conta de poupança na qual os juros são 
continuamente capitalizados numa taxa constante de 5.5% ao ano. A equação y' = ky 
descreve a taxa de crescimento do montante investido, onde k é a taxa de juros. Se $5000 


foram inicialmente investidos, qual é o montante após 3 anos? Resp: $5896, 96 


8. Suponha que a cada mês uma população aumenta na razão k, isto é, no primeiro mês é po, 
no segundo mês kpo, etc. Mostre que a população p(t) satisfaz uma equação diferencial 


do tipo p' = Ap e determine À em função de k. 
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Aula 1.4 


Equações Homogêneas 





Exemplo 1.4.1. 











Xx+y o a E 
(1) F(x, y) = E é homogênea pois 
ctxrty Hrxry) x+y 
Req = aro PO 
e Fl )=1+5=1+u=g(u), onde u = y/x. 
3 2 
(2) F(x, y) = R o é homogênea pois 
Por Plrry) Cray 
Flix, ty) = 58 o > Bm - no - = F(, y) 
Py+rtriy o Plyrxiy) y+xy 
2 
1+(2) 1+u? 
e O a =g(u) onde u=y/x. 


3 
Mt 
x (E) 
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Método de solução 


Equações homogêneas podem ser transformadas em equações separáveis mediante a mu- 


dança de variável natural u = y/x. Desta forma: 


V=ux > í-D=x—+u. 


Substituindo estes dados na equação diferencial, obtemos: 





e pe (u) > pele ()-u > ais 
dx "278 dx” 8 g(u-u x” 
que pode ser resolvido pelo método das equações separáveis. 
dna as : dy 2 
Exemplo 1.4.2. Resolva a equação diferencial ya =y-2. 
Solução: is 
dy dy y— 
2xy- = -2>L=F = 
a UR (x,y) xy 





Da 
A função F(x, y) = — = i E Ea 


Er > |& E 5] é homogênea, logo fazemos a mudança de variável 


u=y/xouy=ãuzx. 


Aplicando a mudança de variável como sugerido acima, chegamos a 





Resolvendo pelo método das equações separáveis, obtemos: 


2udu dx 





+11 x 
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Integrando ambos os lados, chegamos a: 
> 1 
In(u +) =-Inlxj+c=n 5 +c.. 


Aplicando exponencial em ambos os membros, 


P+1=5, onde c=* > 1+(y/32 = À. 
x x 


A solução pode ainda ser reescrita na forma: 


2 


+ 2 CX (x ) + 2 E 
= ou Ei Sis 
y o a Mr 


j c 
Esta solução representa uma família de círculos que passam pela origem e com centro (50) 


sobre o eixo x. 


1.4.1 Translação para Equações Homogêneas 


Uma equação diferencial da forma 


dy E mx+byy+cy 
dx aox + boy + c> 


pode se tornar uma equação diferencial homogênea mediante uma translação (mudança de 


variáveis) conveniente: 


Se os termos constantes cj e c> são ambos zero, então a equação acima já é homogênea. Isso 
fornece a idéia de que transladando x e y por uma determinada constante, a equação resultante 
terá os termos constantes cy e c; ambos iguais a zero. Para transformar esta equação em uma 


equação diferencial homogênea use as substituições: 


x=X+4h 
y=Y+k, 


com h e k escolhidos de tal forma que as duas equações: 


mx+by+c =0, ax+boy+c =0 
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se reduzem às equações 


mX+bY=0, aX+bY=0. 


Portanto escolhemos h e k como soluções do sistema: 


ah+bk+c =0 
ah + bok+ co =0. 
dy x+y+3 


Exemplo 1.4.3. Resolva E E: 


Solução: Resolvendo o sistema 
x+y+3=0 
x->y-5=0, 


temos x = 1, y = —4. Portanto vamos usar a mudança 


x=X+1 
y=Y-4. 


Logo — = no e a equação diferencial se reduz a 


= 


XX y ds 
dx X-Y 1-& 








Y 
Seja U = : ou Y = UX. Então e =U+ oro Usando estas mudanças, temos 





dX dX 
dU 1+U 
ER 


1+U-uUMa-uU 1+uU? 
a ( ) E 





dX es Sei 
f=t dX 
dU=—. 
El X 


Integrando: 


arctg(U) — > Inçl +U?)=InX+Inc 


= 2 arctg(U) = In[cXZ(1 + U?)] 
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> 2 arctg (5) =In [o f + E) 


5 2 arte (2) =In lc (x2 + Y?)| 


> 2arctg (3) =In [c ((x -12+(y+ 9) E 


x—1 


1.4.2 Exercícios 


1. Resolva a equação diferencial dada. 


Pg? 
psd 


(b) ydx+(24X)-x)dy =0 


(c) (12 — y2)dx — 2xydy = 0 


(d) y S 

(e) (1x2 + y)dx + (xº — xy)dy = 0 
dy y 

(£) dx O x+ DT x>0,y>0 





2. Mostre que a equação y' = é homogênea e que possui, para cada ponto (x9, 40) do 


2-y 
plano com yo + xo, uma única solução y(x) satisfazendo y(x9) = Yo. 


3. Resolva a equação diferencial dada. 


dy 2x-y+1 
(a) dx x+3y-2 
dv 8X +34 2 
(b) e 
do =p roy= 
dy  x+3y+4 
MEI 
Ko Ap=54 81 
dy  x+y+3 
(d) dx 2y-x+3 


4. Resolva os seguintes problemas de valor inicial: 
(a) y = y/x+ (2/9) cos (2), vm) = VT. 
(b) yxInx =y,y(3) = In81. 
(o) xy =y+ 4% cos 2(y/x), y(2) = 0. 
() yy =(x-De É, yo) =1. 
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Aula 1.5 


Equações Exatas e Fatores Integrantes 


Seja uma família de curvas f(x, y) = c, onde c é uma constante e f uma função diferenciável. 


Vimos do Cálculo que a diferencial total de f é: 


mB as O 
df = Eid + ou 


Disto segue que se f(x, y) = c, então df = 0. 


Exemplo 1.5.1. Se f(x,y) = x2yº = c, então a diferencial total desta expressão é: 
df = 2xy'dx + 3x22yºdy = 0. 


Reescrevendo a equação acima e dividindo tudo por dx, obtemos: 





dy Dye 
dx 3x2y” 
Isto é, considerando f uma função diferenciável, a uma família de curvas f(x,y) = c sempre 
9 o 
corresponderá uma equação diferencial do tipo Lar + Srdy = 0. 
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Observação: Segue da definição que se M(x, y)dx + N(x, y)dy = O é exata, então por (1.5.1) e 


(1.5.2), existe uma função F(x, y) tal que 


oF oF 
a Mt(x, y), Po N(x, y) 


e F(x,y) = cé a solução geral da equação ( 1.5.1). 


Exemplo 1.5.2. Seja a equação diferencial 


(9x? + 2yº + 2) dx + (4xy + 12º) dy = 0. 
e — me” e — ee 


M(x,y) N(xy) 
4 ns 3 - 
Ela é exata. Note que se F(x, y) = 3x2 + 2xyº + 4º + 2x, então 


OF 


er = (92º + 2y +2)=M(x,y) e e (4xy + 124º) = N(x,y) 


ox 
Logo 3% + 2xyº + 4yº + 2x = c é a solução geral da equação diferencial. 
fed y 8 q 


Método de solução 


Antes de descrevermos o método de solução, precisamos analisar duas questões importantes 


com relação às equações exatas. 


1. Existe uma maneira de testar se uma equação é exata ou não? 


2. Se uma equação é exata, como poderemos achar sua solução geral, isto é, como poderemos 
encontrar uma função F(x, y) tal que F satisfaz a equação exata e F(x,y) = c é a solução 


geral? 
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3. Se uma equação não é exata, é possível modificá-la para torná-la exata? 


1.5.1 Determinando se uma equação é exata 
Se F(x, y) = c é uma solução da equação diferencial exata M(x, y)dx + N(x, y)dy = 0, então 


OF 


a M(x,y) e Cl N(x, 4) (1.5.3) 


oy 
Suponha que Me N são contínuas e com derivadas parciais contínuas num domínio retan- 
gularR: a<x<B, y<y<oódo plano. Calculando a derivada parcial de M em relação a y e 


de N em relação a x, obtemos: 





Pr OM 
dyvox Oy 
92F ON 
dxdy Ox 
Pela hipótese de continuidade das derivadas parciais M, e Nx, segue que as segundas derivadas 
Gi ia gi cai a poa RD CIO O Aa 
parciais mistas são iguais EE RR isto é, 
oM ON 
o ar (1.5.4) 


Resumimos a ideia acima no seguinte teorema. 


oN dM : 
Teorema 1.5.1. Sejam as funções M,N, x e Em contínuas num domínio retangular 


R:a<x<bB, y<y<ódoplano. Então a equação diferencial M(x, y)dx + N(x, y)dy = O é exata 


em R se, e somente se, 


om ON 
dy dx 





Exemplo 1.5.3. Verifique que a equação diferencial cos (x+ y)dx+(3yº +2y+ cos (x+y))dy = 0 


é exata. 


Solução: Note que M(x,y) = cos (x+y) e N(x,y) = (3yº + 2y+ cos (x + y)). Se a equação é 
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É oM ON 
exata, então devemos ter — = — De fato, 
dy dx 


Jy sen y 
Ba = — sen (X + y). 


Portanto, a equação é exata. 


1.5.2 Determinando a solução geral 


Se a equação diferencial ( 1.5.1) é exata, queremos determinar uma função F(x, y) tal que 
dF(x, y) = M(x, y)dx + N(x, y)dy, 


com 


oF 
sia M(x, y) 
(1.5.5) 
OF 
om = N(x,y) 
Integrando a primeira equação de ( 1.5.5) em relação a x, obtemos: 
OF 
[5a = fra, ydx 
es F(x,y) = fa, y)dx + k(y), (1.5.6) 


onde k(y) é uma constante de integração quando se integra em relação a x. Para determinar 


k(y), derivamos a função F em relação a 1 e igualamos a N, como na equação (1.5.5): 


OF 9 ; 
Em = Er [ [ma pla) +K(y) 
es k(y) = N(x,y)- s/a. pj) (1.5.7) 


Esta última equação determina a função k(y). 


Exemplo 1.5.4. Resolva a equação diferencial cos (x + ydx + (3yº + 2y + cos (x + y)dy = 0. 


MATB97 - Equações Diferenciais 


43 


Solução: Já vimos no exemplo acima que esta equação diferencial é exata com 
M(x,y) = cos(x+y) e N(x,y) = (3º +2y+ cos (x + y)). 


OF 
Agora precisamos determinar sua solução, isto é, uma função F(x, y) tal que — = cos (x+y) e 


Ox 
OF 
eo = N(x, y). Integrando — = cos (x + y) em relação a x, obtemos: 


Oy Ox 
F(x,y) = cos (x+y)dx+k(y) o F(x,y)= sen(x+y)+k(y). 
Derivando a função F obtida em relação a y e igualando a N(x, y), obtemos: 
3 = cos(x+y)+k(W=3W+2y+cos(x+y) é kKy=3+y o Ky=yp+y+a, 


onde cj é uma constante real. Substituindo k(y) na expressão de F acima, encontramos 


F(x,y)= sen(x+y)+y +y +c1 e, portanto, a solução geral é: 
sen(x+)+y+y+a=o es sen(xty)+yp+y =, 


ondec=c,—41. 


Exemplo 1.5.5. Acheo valor deb para o qual a equação diferencial (xy? +bx? y)dx+(x+y)x2dy = O 


é exata e resolva a equação com este valor de b. 


Solução: 
Note que 
Mx, y) =x +by e N(x,y)=(x+y. 
Para a equação ser exata, precisamos ter a = E isto é, 
dy dx 


2xy+ br =3x2+2xy O b=3. 


Assim, vamos encontrar a solução de (xy? + 3x2ydx + (x + yx2dy = 0, isto é, encontrar uma 


OF 
função F(x, y) tal que - = M(x,y) e 5 = N(x,y). Integrando o xy? + 8x2y em relação a x, 
obtemos: 

2d 


F(x, y) = fe +32 y)dx = — +êy + k(y). 
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OF 
Derivando esta última expressão em relação a y e observando que — = N(x, y), temos: 


dy 


= dy +% +K'(y) = + yx? > ky=c. 


Substituindo k(y) na expressão de F, a solução geral será: 


1.5.3 Fator integrante: transformando uma equação diferencial em 


exata 


Exemplo 1.5.6. A equação diferencial -ydx + xdy = O não é exata. De fato, 


oM 
Mx, y)=-y > au 
N(xy)=x > q. 


o 
Como a + a? equação não é exata. Porém, note que se a equação diferencial for multipli- 
cada por 1/x2, ela se torna 


a + 7iy =0, 


que é uma equação exata: 


Mer > dy 1/x? 
N(xy)=1/x > Day e 


e como —— = —=—, esta nova equação diferencial é exata. 
y x 


Este último exemplo nos mostra que, em geral, se a equação diferencial 
M(x, y)dx + N(x, y)dy = O (1.5.8) 


não é exata, podemos multiplicá-la por uma função apropriada (x,y), de modo que a nova 
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equação diferencial 


M(x, ul, y)dx + N(x, y) lx, ydy = O (1.5.9) 


seja exata. 
A função u(x, y) é chamada de fator integrante da equação ( 1.5.8). A equação ( 1.5.9) pode 
então ser resolvida pelo método da equações exatas e suas soluções também serão soluções da 


equação original (1.5.8). 


a Ia 


Exemplo 1.5.7. Verifique que u(x, y) = — é um fator integrante para a equação diferencial 


ydx+2xInxdy=0 
e então encontre sua solução. 


Solução: Seja M(x,y) =y e N(x,y) = 2xInx. Então 


EM ue e ON minas 
Oy Ox 


oM | ON 
Como — + — a equação dada não é exata. Multiplicando a equação por - temos a nova 


dy dx 
equação 
y 
E dx +2ylnx dy = 0 
Agora se M(x, y) = e N(x,y) = 2yln x. Então 
0M 2y , ON (24 
du Cdx 
oM ON ai Hs ; 
Como or = x à Nova equação é exata e portanto u(x, y) = ; é um fator integrante da equação 
diferencial original. Como a nova equação é exata, existe uma função F(x, y) tal que 
EE 
dx x 
(1.5.10) 
5 =2ylnx 
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OF ; 
Integrando —— com respeito a x, temos: 


Ox 


F(x,y) = y? In x + h(y) 


Derivando F(x, y) em relação a y e igualando a N(x, y), temos: 


OF ; 
Em =2ylnx+h'(y) 


> h(y=0=>h(y=c c= constante. 


Portanto, substituindo h(y) na expressão de F, a solução geral será: 
ping= é 


Determinando fatores integrantes 


Seja M(x, y)dx + N(x, y)dy = O uma equação diferencial não exata com fator integrante (x, y). 


Então 


M(x, y)u(x, dx + N(x, Wulx, y)dy = 0 (1.5.11) 
seja exata. Neste caso deveremos ter 


oMu  ONu ou oM du ON 
a CR o a PONTE hm 
ou 


2 a 
É E [5 a ] (1.5.12) 


=—M — =-N = DES o 

Ou seja, se pudermos encontrar a função (x,y) satisfazendo a equação (1.5.12) então a 
equação ( 1.5.11) será exata. Infelizmente a equação ( 1.5.12), que determina o fator integrante é 
muito difícil de ser resolvida. Desta forma, procuramos fatores integrantes nos seguintes casos 


especiais: 
Primeiro caso: | é uma função somente em x, isto é u = (x) 


o 
Se u = u(x), então Cefei equação ( 1.5.12) se reduz a: 


dy 


d 
En-(S | 


dx Oy Ox 


dy dx E 
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ou 


9M — ON 
du E o 
E = RR ul). (1.5.13) 
a — 


função só em «x 


oM — 9N 


o 
Mas a equação ( 1.5.13) não tem sentido, a menos que a expressão RE seja uma função 


somente em x. Isto é, 
(2 2N) 


Oy ox 


NT = (a) 


para alguma função h de uma variável, e a equação ( 1.5.13) se reduz a 


du 


die at 


que é uma equação de variáveis separáveis. Separando as variáveis temos que 


e integrando temos que 


ny = [no dy => q = elo) dx, 
Concluímos então que a equação admite o fator integrante u(x) : 
u(x) = ele) dr, (1.5.14) 


Segundo caso: | é uma função somente em y, isto é u = u(y) 


9 
Se u = u(y), então = 0 ea equação ( 1.5.12) se reduz a: 


d 
Em = (5 e) 


dy dy Ox 
Rs Em a) 
du O dy dr 
a ME u(y). (1.5.15) 
O — 


função só em y 
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oM — ON 


o o 
Mas a equação ( 1.5.15) não tem sentido, a menos que a expressão E seja uma função 


somente em y. Isto é, 


Em E a) 
Oy Ox 


RR 


para alguma função f de uma variável e a equação ( 1.5.15) se reduz a 


d 
ay Tt fO) 


que é uma equação de variáveis separáveis. Separando as variáveis temos que 


e integrando temos que 


ng = - [rm dy p= Say, 


Concluímos então que a equação admite o fator integrante u(y): 


u(y) = e! fúnay, (1.5.16) 


Terceiro caso: u é uma função de xy, isto é, (x,y) = u(xy) 


9 o 
Se u(x, y) = u(t), com t = xy, então — =y-u(t) e 5 =x-u'(t). A equação (1.5.12) então 
se reduz a: 
inte [o en 
x: (OM y- WON = [5 =) () 
e (2 2N) 
du O Nov dk 
O — 


função só em + 


oM — ON 
QM ON 
Mas a equação ( 1.5.17) não tem sentido, a menos que a expressão MON seja uma função 


somente em t. Isto é, 


(4-8) 
MON é 80 
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para alguma função g de uma variável, e a equação ( 1.5.17) se reduz a 


du 
dE uso) 


que é uma equação de variáveis separáveis. Separando as variáveis temos que 


e integrando temos que 


Inu =— EO dt = u(t) = e JS st) dt 


Conclui-se então que a equação admite o fator integrante u(xy): 
u(t)=e Íst)dt, onde t=xy. (1.5.18) 


Resumimos a teoria dos fatores integrantes no seguinte Teorema. 


Teorema 1.5.2. Seja M(x, y) dx + N(x,y) dy = O uma equação diferencial não exata. 


oM — ON 
dy dx 
e Se [= 5) = h(x) ( é uma função somente em x), então a equação admitirá o fator integrante 
u(x) = eJ o) dx 
(84-83) 
o Se Re = f(y) (é uma função somente em y), então a equação admitirá o fator integrante 
u(y) = e JS fe) dy 
(ut -28) 
e Se E = g(t), com t = xy, então a equação admitirá o fator integrante 


ut) = e Su) di 


Exemplo 1.5.8. Resolva o PVI 


2 sen (y?)dx + xy cos (y')dy = 0, yv2) = vn/2. 
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Solução: Vamos testar se a equação diferencial é exata: 


M(x, y) = 2 sen (y?) > 5 = 4y cos (7) 
N(x,y) =xy cos(y?) > = =ycos (y?) 
Como ua f eu a equação não é exata 
o om Srta 
Fator integrante: Temos que 
oM ON E- 2 > o 2 
E 5) = arcos (1) y cos (yº) = 3y cos (yº). 
Assim 
(g- 28) 
dy dx) 3 
N RE 


Pelo Teorema 1.5.2, a equação admite o fator integrante 
u(x) = eJii=dno u(x) = xº. 


Multiplicando este fator integrante pela equação, obtemos 2x? sen (y?)dx + x*y cos (y?)dy = 0. 


Vamos testar se a nova equação diferencial é exata: 


M(x,y) = 2% sen (y)) > Em = 422 y cos (y?) 
4 ON 3 2 
N(x, y) = xy cos (y) > a y cos (yº) 
oM ON e é : = 
Como E E a equação é exata. A seguir, encontraremos a solução geral F(x, y) = c, onde 
e =Me cb = N. Integrando as = M em relação a x, obtemos: 
Ox Oy Ox 


F(x,y) = [ 2% sen (y?)dx + k(y) = = sen (y?) + k(y). 


De 3 = N, temos: x!y cos (y) + K(y) = x*y cos (7) => Ky) = 1. Logo, a solução geral 


1 
será F(x, y) = c, isto é, 5“ 


sen (1?) = c. Substituindo a condição inicial, y = 7/2 parax = 2: 


1 
516 sen (r/2)=c 5 c=8. 
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A 


1 
Portanto, a solução do PVI é portanto 2” sen (W) =8. 


Exemplo 1.5.9. Resolva a equação 
(x2 + yº) dx + (xº +3xyº + 2xy) dy = 0 


Solução: Vamos testar se a equação diferencial é exata: 


Mxy)=2+y > o =2y 


N(x,y) = +3xy2 +2xy > map +dy 


oM ON ipa coa 
Como —— * —, a equação não é exata. 


dy dx 


Fator integrante: Temos que 


oM ON E e 2 2 ES a 2 2 
E e) = 24 (3xº +3yº +2y) = —3(xº + yº). 
Assim 
oM ON 
(4-5), 
M =—3. 


Pelo Teorema 1.5.2, a equação admite o fator integrante 
u(y) = e SB =, 
Multiplicando a equação pelo este fator integrante, obtemos 
(x? + yº) dx + Y(xê + 3xy? + 2xy) dy = 0. 


Vamos testar se a nova equação diferencial é exata: 


oM 


M(x,y) = Ex? + y?) > 27 = 34x? + y?) + 2ye4 = (3x? + 3y? +2y) 
34/13 2 ON ay 2 
N(x, y) =x +3xy +2xy) > a Te + 3yº + 21) 
oM ON o a E 
Como E = E a equação é exata. A seguir, encontraremos a solução geral F(x, y) = c, onde 
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top 


dE M OF 


OF 
e — = N. Integrando — = M em relação a x, obtemos: 
Ox Oy 


Ox 


F(x,y) = fere + y?) dx + k(y) = ey [5 + nº + k(y). 


De E N, danada EN (x? ES 3xy? dE 2xy) + K(y) = E (é E 3xy? + 2xy) >> K(y) =0 => 


3 
k(y) = c1. Logo, a solução geral será F(x, y) = c, isto é, ey [5 + ui =coux + 3xy? = ce. 


Exemplo 1.5.10. Resolva a equação 
É 
[ox+ axe (Eu !) dy=0 
b y 


X 


Solução: Vamos testar se a equação diferencial é exata: 


6 oM 6 
Mw,yj=3k+—- > D=-S 
á y dy y 
e NO ON 2x 3y 
Ne a E 
C oM ON ão não é xat 
omo y E ESA AÇE O ao e exata. 


Fator integrante: Temos que 


oM ON - a) E = =] —6x? — 2x3y + 3y? 


dy dx y v 2 2y? 
E 2 3 3 
2 3) 62+28y-3 
um = qn = (314 6) y(5 + Tt) OE 
y vox xy 
Assim 


oM ON 

(54 2) 1 1 
ND =D-- = g(t). 
xM — yN xy t 


Pelo Teorema 1.5.2, a equação admite o fator integrante 
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Multiplicando a equação pelo este fator integrante yu = xy, obtemos 
(3x2y + 6x) dx + (1º + 3y7) dy = 0. 


Vamos testar se a nova equação diferencial é exata: 


M(x,y) =3xy+6x > GM anão 
dy 
ON 
-.3 2 = 242 
Ne, yW)=02+3y > DK 3x 
oM ON Eu E : a 
Como Br E a equação é exata. A seguir, encontraremos a solução geral F(x, y) = c, onde 
OF oF OF 
— =Me — =. Integrando — = M em relação a x, obtemos: 
Ox Oy Ox 


F(x,y) = [xy + 6x) dx + k(y) = 22y + 3x2 + k(y). 


De 5 = N, temos: é +K'(y) — mw +3y? Es K(y) = 34º sá; k(y) E P Per Logo, a solução 


geral será F(x,y) = cisto é, x22y +32 + y' =c. 


1.5.4 Exercícios 
1. Verifique se cada uma das equações a seguir é exata. E, caso afirmativo, resolva-a. 


(a) (2x +3) + (2y —-2)y' = 0. 
Resp: x2+3x+yº -2y=c. 

(b) (3x2 — 2xy + 2)dx + (6º — x2 + 3)dy = 0. 
Resp: 12º —- x2y+2x +22 +3y=c. 
dy — ax+by 

(ao bx+cy 
Resp: ax? + 2bxy + cy? = k. 

dy — ax-—by 

pm bx— cy 

Resp: Não é exata. 





(e) (e sen y — 2y sen x)dx + (e cos y + 2 cos x)dy = 0. 
Resp: e seny+2ycos x=ce também y = 0. 
(9) (ye cos 2x — 2º” sen 2x + 2x)dx + (xe” cos 2x — 3)dy = 0. 


Resp: e! cos2x+xº -3y=c. 
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(g) (xIny+xydx+(yinx+xy)dy=0 x>0,y>0. 


Resp: Não é exata. 
2 2 dy 
(h) xy + y+(x +) =0 


dy 


(i) € sen y+2x+ (eº cos y +24) — =0 


d 
0) ML + xy) + rea =) 


(1) (2x cos y—- eJdx — x2 sen ydy =0 
(m) (3x2 y? + sen x)dy + (2x? + y cos x)dx = 0 
2. Resolva os seguintes problemas de valor inicial 
(a) (2x — y)dx + (Qy — x)dy = 0, y(1) = 8. 
Resp: xº -xy+yº =7. 
(b) (9x2 + y — 1)dx — (4y — x)dy = 0, y(1) = 0. 
Resp: 32 +xy-x-2yº =2. 


3. Mostre que qualquer equação separável, 
M(x) + N(y)y' = 0, 


é também exata. 


4. Mostre que as equações a seguir não são exatas, mas se tornam exatas quando multipli- 


cadas pelo fator integrante dado. A seguir, resolva cada equação. 


(a) Pp rx +y)y' =0, u(x,y) = 1/xyº. 
Resp: x? +2Inly- y? = cetambém y =0. 

(b) [= — 2e “ sen ) dx + E dv=0, uism=yve. 
Resp: é seny+2ycosx=c. 

(c) ydx + (2x — ye)dy = 0, (x,y) =y. 
Resp: xy? — (yº —-2y+7)e! =. 

(d) (x+2) senydx+x cos ydy=0, u(x,y) = xe. 


Resp: x2e” sen y = 0. 
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(e) (2 +y -xdx-ydy=0; ulxy=(C+y)! 
5. Para as equações diferenciais abaixo, ache um fator integrante e resolva a equação. 


(a) (322y + 2xy + y)dx + (1º + x2)dy = 0. 
Resp: u(x) = e, (3x2y + y Je” = 
b) y=+y-1. 
Resp: u() =, y=c+1+e2. 
(c) dx + (x/y — sen y)dy = 0. 
Resp: u(y) = y xy+ycosy-— seny=c. 
(d) ydx + (2xy — e 2Ndy = 0. 
Resp: u(y) = e /y, xe! — In |y| = ce também y = 0. 
(e) Cdx + a cot y + 2ycsc y)dy = 0. 
Resp: u(y) = senye seny+y =c. 
o (ax+ DA TER “lã =0. 
vox 
Resp: u(x,y) = xy Xy +32 +y' =c. 
(g) x? +2y-— = =0 


(h) xe" + xlny + e aa JE -o 
y+y y AE 


(1) ydx + (2xy- e 2YN)dy =0 


4º 3 3x 
(j (+ ae (Sea Ja =0 
) o Rê y|dy 
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Aula 1.6 


Equações Lineares de primeira ordem 





Ê 
Exemplo 1.6.1. A equação diferencial y' + ida lélinear comp(x) =1/2eg(x)=1. 


Exemplo 1.6.2. A equação diferencial y +5y = sen x não é linear. 


1.6.1 Método de solução 


A equação diferencial linear de primeira ordem 


d 
E +peyy=40) (1.62) 
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pode ser resolvida com ajuda de um fator integrante. Reescrevendo ( 1.6.2) na forma 


M(x, y) dx + N(x,y) dy = O 


isto é 
(p(xc)y — qg(x)) dx + dy = 0 (1.6.3) 
temos que 
oM 
M(x, y) = pO)y — qu) = air pl) 
N(g,y)=1 > D=0 


oM ON 
Como — + —, a equação ( 1.6.3) não é exata. 


dy dx 


Fator integrante: Temos que 


(e -) nto 


dy Ox 
Assim 
oM ON 
e) 
EN RAS) 


Pelo Teorema 1.5.2, a equação ( 1.6.3) admite o fator integrante 
ua) = elo dx, 
Multiplicando a equação ( 1.6.2) por este fator integrante, obtemos 
x 


ed pt) dx, Ea pl)y . eJ re) dx — g(x) E el P(9 dx (1.6.4) 


O fato importante da equação ( 1.6.4) é que o lado esquerdo pode ser escrito como 


d (y- el vt ca) 


d 
[px da, LY cof) dx — A 
e EE + p(x)y-e EE 


Portanto a equação ( 1.6.4) ficará 


d (y- el PO 2%) = q(a9 - el rt ax 


dx 
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Integrando temos 


yreJoeide — [qq velo de duda 


Portanto a solução da equação ( 1.6.2) é 


y= e S Pto) dx TES cel qr+c (1.6.5) 


Teorema 1.6.1. À equação diferencial linear de primeira ordem 


y + pl)y = q(x) 


admite um fator integrante u(x) = e! PC) dx e q sua solução geral é 


y = e Sl) dx | faoo «el re) dx dx+c 





Observação. Não é necessário memorizar esta fórmula. Para resolver a equação diferencial linear, 


y + plOdy = q() 
(1) calcule o fator integrante u(x) = el PG) dx, 


(2) multiplique a equação pelo fator integrante e isto transformará o lado esquerdo da equação em 
d 
(y celrto dm; 


dx 


(3) integre para achar a solução geral. 
Exemplo 1.6.3. Resolva y' + - =3c0s2x, x>0. 


Solução: A equação diferencial é linear. 


: 1 
Fator integrante: eli = ent = 


Multiplicando a equação pelo fator integrante temos: 


xy +y=3x cos 2x. 
Nem 


É (xy) 
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Portanto, 


d 
ax tv) = 3x cos 21. 


Integrando, temos: 


Xy = 3 fxcos 2x dx 


(integração por partes) 
1 1 
ay = 3 [grsen Ze 5] sen rca o 
1 1 
5 axy = 3 [5x sen 2x+ Z cos 2x) +c 


Portanto, a solução geral é: 


SE DO cos 2x c 
Po 4 





X 


Exemplo 1.6.4. Ache a solução do problema de valor inicial y +2y = e“ ey(0) = 3/4. 


Solução: A equação diferencial é linear. 
O fator integrante será: 


u(s) E e) 2dx Es e, 


Multiplicando a equação pelo fator integrante temos 


E y +24 =. 
Ven mem” 


E (ey) 


Portanto, 


Integrando, temos: 


Portanto, a solução geral é: 
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Substituindo a condição inicial: 3/4 = c+1. Logo c = —1/4e a solução do PVI é 
y=-0.25 * +e*, 
Exemplo 1.6.5. Resolva a equação diferencial xy' + (x +1)y =x. 


Solução: A equação diferencial não está na forma de uma equação linear. Precisamos dividir 








+1 +1 
tudo por x, obtendo a equação equivalente y' + E p= na forma linear com p(x) = E 
Pp q q Y EA p E 
g(x) = 1. O fator integrante será: 
u(s) = eJ(I+i)dx = gvtinx = xe". 
Multiplicando a equação pelo fator integrante temos: 
xety' + (x + 1)e'y = xe 
Vecemcanenssamm meme” 
de (ney) 
Portanto, 
que 4) Ed cdi 
Integrando, temos: 
ME = HR xe” dx 
(integração por partes) 
>xey = xe -fe dx 
>»xedy = x -e+c. 
A solução geral é: 
E ge da 
PEC pe 
que é a solução geral da equação diferencial linear. 
Exemplo 1.6.6. Seja y = yi(x) uma solução da equação 
y +p()y =0 (1.6.6) 
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e seja y> = y>(x) uma solução de 


y + pO)y = q(9). (1.6.7) 
Mostre que y = Wj + y> também é solução da equação ( 1.6.7). 


Solução: Para mostrar que a função y = y + y> é solução da equação ( 1.6.7), devemos mostrar 


que ela satisfaz a equação diferencial. De fato, 


(1 + 97)" + pl +42) = q(o) 
W HP + +pldy = quo) 
O+y+plyp = q(x) 

qo) = q(9). 


1.6.2 Exercícios 


1. Ache a soluçao geral das seguintes equações diferenciais lineares de primeira ordem: 


(a) y +3y=t+e?. Resp: y=ce+(t/3)-1/9+e 2 
(b) y +y=te+1.Resp:y=ce!+1+te'/2. 

(o) y —2y = 3€. Resp: y = ce! — 3. 

(d) y +2ty = 2te É. Resp: y = Pef + cof. 

(e) 2y +y=3t. Resp: y=ce!2+3t- 6. 

(f) y +y=5sen2t. Resp: y=ce! + sen2t-2 cos 2t. 
(9) W — (tg Dy = cost. 


du xe +u 
(h) = 


dx x 
(D) y 





2dy + ydx = xdy 
D y -4y=2x- 4% 
(k) xdy + (y — sen x)dx = 0 
2. Resolva os problemas de valor inicial abaixo, 
(a) y — y=2te!, y(0) = 1. Resp: y = 3€ + 2(t — 1)e2. 


(b) ty +2y=P-t+1, y(1)=1/2,t>0. Resp: (3! — 4Ê + 62 + D)/12Ê. 
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(O) y —2y= 2, v(0) = 2. Resp: y = (t+ 2). 
(d) By +42y=e!,y(-1)=0. Resp: y= (t+ De ff. 


(e) y —(tgDy= cost, yv(0) = 1. 


dy 
(8 o Vig = sec t 
y(O) = O 
E a 
(g) dx x 
yle) = O 


3. Considere o PVI 
o Jd 


Ache as coordenadas do ponto de mínimo local da solução. Resp: (t,y) = (-2In - —4ln 5) 


4. Considere o PVI 


Ud 1 
y Faves, y(O) = yo. 


Ache o valor de yo para o qual a solução toca, mas não cruza, o eixo Ot. Resp: yo = 


—1.642876. 
5. Ache o valor de yo para o qual, a solução do PVI 
y—-y=1+3sent, y(O)=yo 
permanece finito quando t — oo. Resp: —5/2. 
6. Mostre que se a e À são constantes positivas e b é um número real qualquer, então toda 


solução da equação 


y +ay = bet! 


tem a propriedade que y — 0 quando t — oo. Sugestão: considere os casos a = À e 


a £ À separadamente. 


7. (Equação de Bernoulli:) Seja a equação 


y = fO)y + go)y, (1.6.8) 
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onde « £ 0,a + 1 é uma constante e f e q são funções conhecidas. Faça a mudança de 


variável u = y!”*. Mostre que esta substituição reduz a equação acima à equação linear: 
uv =(1-0)f(x)u + (1 — a)g(x). 


8. Use o exercício ( 1.6.8) para resolver as seguintes equações: 


4y 2.4 4 E 2 
(a) Ra y. Resp: y=cix*+cx Inx + (nal. 


fo Ds maço: Y ; — ca 
(b) y' = —xy” + Resp: y RO 
(ey + =y =y? cosx,y(l)=1. 

A 1 
(d) y + 2y= uy, >) =—2 


9. Seja y = yi(x) uma solução da equação 
y +pl)y=0 (1.6.9) 
e seja y = y>(x) uma solução de 
y + plJy = g(x). (1.6.10) 
Mostre que y = y1(x) + y>(x) também é solução da equação ( 1.6.10). 


10. Mostre que a equação diferencial y' + p(x)y = g(x)yln y pode ser resolvida usando a 


mudança de variável z = Iny. Aplique este gmétodo para resolver a equação xy' = 


cx+2x2 


22y+ ylny. Resp: y=e 


11. Mostre que a equação ( cos y)y' +2x sen y = —2x pode ser transformada em uma equação 


linear com a mudança de variáveis z = sen y. 


12. (Variação de Parâmetros): Considere o seguinte método para resolver uma equação 


diferencial linear geral de primeira ordem: 
y + p(Dy = g(). (1.6.11) 
(a) Se g(t) = 0, mostre que a solução é: 


y = Ae Jr (1.6.12) 
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(b) Se g(t) 0, suponha que a solução é da forma 
y=A(ge Srbat, (1.6.13) 


onde A é agora uma função de t. Substituindo y e y' na equação diferencial ( 1.6.11), 


mostre que A(t) deve satisfazer a condição 
A'(b) = g(Del PO, (1.6.14) 


(c) Ache A(t) a partir de (1.6.14). Substitua o valor encontrado na equação ( 1.6.13) e 
obtenha a solução da equação diferencial. Este método é conhecido por variação de 


parâmetros. 


13. Use o método de variação de parâmetros para resolver a equação diferencial y —2y = Pe, 


t>0. Resp: y(t) = cl! + Pela. 


14. (Equação de Ricatti): A equação 


v + pl)y + go) = f(x) (1.6.15) 
é conhecida como equação de Ricatti. 


(a) Mostre que se y1(x) e y>(x) são soluções da equação ( 1.6.15), então a função z(x) = 


y>(x) — yi(x) é solução da equação de Bernoulli 
Z+(p+2ygz-q7 =0. (1.6.16) 
(b) Sabendo que y(x) = x é uma solução de Ricatti 
y + êy — Re =1, 


determine as demais soluções. 
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Aula 1.7 


Aplicações das Equações Diferenciais 


de primeira ordem 


Modelagem é o processo de escrever uma equação diferencial para descrever uma situação 
física. 

Muitos fenômenos em ciência e engenharia envolvem relações entre quantidades que variam 
com o tempo. E como taxas de variações são representados matematicamente por derivadas, 
tais fenômenos podem ser modelados por equações que envolvem derivadas. Nesta aula vamos 


analisar a modelagem de alguns destes fenômenos: 


1.7.1 Problemas de variação de temperatura 


A Lei do resfriamento de Newton afirma que a “taxa de variação da temperatura de um corpo 
é diretamente proporcional à diferença entre a temperatura do corpo e a temperatura do meio 
ambiente.” 


Modele este fenômeno como uma equação diferencial de primeira ordem e classifique-o. 
e Seja T(t) a temperatura do corpo no tempo te T,, a temperatura do ambiente. 
O ds .dT 
e Ataxa de variação da temperatura é ar 
e Então o modelo matemático para a lei de resfriamento de Newton é: 


dT 
dt = KT > Tm), (1.7.1) 


onde k é um constante de proporcionalidade. 
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e A equação (1.7.1) pode ser resolvida como uma equação diferencial linear , pois podemos 


reescrevê-la como 
dT 
— =kT==hT,. 
dt E 


ou como uma equação diferencial de variável separável, pois podemos separar os variá- 
veis: 
dT 


T-T, =kdt. 





1.7.1.1 Exercícios 


Use a Lei do Resfriamento para responder às seguintes questões: 


1. Um corpo com temperatura desconhecida é colocado em um refrigerador mantido à 
temperatura constante de 0ºF. Se após 20 minutos a temperatura do corpo é 40ºF e após 


40 minutos é 20ºF, determine a temperatura inicial do corpo. 


2. Um corpo à temperatura de 50ºC é colocado em um forno cuja temperatura á mantida 
à 150ºC. Se após 10 minutos a temperatura do corpo aumentou em 25ºC, determine o 


tempo necessário para o corpo atingir os 100ºC. 


3. Ocorpo de uma vitima de assassinato foi descoberto. O perito da policia chegou à 1:00h da 
madrugada e, imediatamente, tomou a temperatura do cadáver, que era de 34,8ºC. Uma 
hora mais tarde ele tomou novamente a temperatura e encontrou 34,1ºC. A temperatura 
do quarto onde se encontrava a vítima era constante a 20ºC. Estime a hora em que se deu 
a morte, admitindo que a temperatura normal de uma pessoa viva é 36,5ºC. 


Resp: t = 2,24 horas. 


4. Um termômetro é retirado de uma sala e colocado do lado de fora em que a temperatura 
é de 10ºC. Depois de 1 minuto a leitura do termômetro é de 15ºC e após 2 minutos 12ºC. 
Determine qual a temperatura da sala onde se encontrava o termômetro inicialmente. 


Resp: 40ºC. 


1.7.2 Problemas de mistura de líquidos 


Há um certo número de exemplos comuns que podem ser classificados como “problemas 


de misturas” cuja configuração é a seguinte: 
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e Temos um reservatório (piscina, lago, oceano) de líquido (água, gás) que tem alguma 


substância (poluição, sal) dissolvido nele. 


e O reservatório começa com um volume inicial Vo e há uma quantidade inicial de subs- 


tância yo no reservatório. 


e Temos uma determinada quantidade de líquido E(t) entrando no reservatório com um 
certo concentração k(t) da substância e uma determinada quantidade de líquido saindo 


s(). 


e Partimos do princípio de que a substância é uniformemente e perfeitamente misturados 
no reservatório, podemos perguntar a quantidade y(t) de substância que permanece no 


reservatório após tempo t. 


e Evidentemente, a quantidade de líquido entrando ou saindo podem ser constantes (ou 
seja, não dependem de tempo), e de igual modo a concentração do líquido entrando 


também poderia ser uma constante. 


Modele este fenômeno como uma equação diferencial de primeira ordem e classifique-o. 


e Seja V(t) o volume total do reservatório em tempo t, com volume inicial Vo 


e Seja y(t) a quantidade total da substância no reservatório no tempo t, com quantidade 


inicial yo 


e Como a concentração da substância entrando no reservatório é k(t) e taxa de entrada do 


líquido no reservatório é E(t), então 


a taxa de entrada da substância taxa de entrada concentração da substância 


no reservatório no tempo t do líquido no liquido na entrada 


= E()-K(t) 


: aê esa) 
e Do mesmo modo, a concentração da substância no reservatório é >— e como a taxa de 


V(t) 
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saída do liquido no reservatório é S(t), então 


a taxa de saída da substância taxa de saída concentração da substância 
no reservatório no tempo t do líquido no liquido na saída 
Ê 
= Ss(t)- ve) 
V(t) 


e Portanto a taxa segundo a qual y(t) está variando com o tempo será a taxa líquida: 


dy a taxa de entrada da substância a taxa de saída da substância 
dt no reservatório no tempo t no reservatório no tempo t 
Pp p 
y(t) 
= E(D)-kb— S(t). — 
O-KD-SO- 7 


e Reescrevendo a equação na forma padrão temos: 


dy Do. 
dE vo O ER) 


que é uma equação diferencial linear. 


Exemplo 1.7.1. Um tanque com capacidade de 400 L contém inicialmente 100 L de água 
salgada com 50 kg de sal. Água salgada contendo 1 kg de sal por litro entra no tanque a uma 
taxa de 5 L/min, e a água misturada no tanque flui para fora a uma taxa de 3 L/min. Quanto 


sal conterá o tanque quando estiver cheio de água salgada? 


Solução: Seja y(t) a quantidade se sal no tanque no tempo t. 





a taxa de entrada de sal taxa de entrada concentração de sal 
no tanque no tempo t de água salgada no fluxo de entrada 
5L lkg É 
E (1) [5Ê) Ro 


Como água salgada entra no tanque a uma taxa de 5 L/min e sai a uma taxa de 3 L/min, em 


tempo t teremos um volume total de (100 + 2t)L de água no tanque. 
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a taxa de saída de sal taxa de saída concentração de sal 


no tanque no tempo t de água salgada no fluxo de saída 





E E Di rã [= 


Imin/ (100 +20)” Re 


Portanto a taxa segundo a qual y(t) está variando com o tempo será a taxa líquida: 





dy a taxa de entrada de sal a taxa de saída de sal 
dt no tanque no tempo t no tanque no tempo t 
z Su 
n 100 + 2t 
ou 
dy 3 
+1———— .y= ki 
di nos O o 
EI = e) mex dt! = e3 In(000+) — (100 + 2832 
Logo 


- | y(100 + 282] = 5(100 + 242 


Integrando temos que 


v(100 + 2t)º/2 


5 [ (100 + 232 dt 

5(100 + 24)/2 
5/22 

(100 + 24)/2 + C 


Portanto, 


y = 100 + 2t + C(100 + 28) */2. 


Quando t = 0, y = 50 => 50 = 100 + €(100)*/2 = C = -50000. 
Logo, 
y = 100 + 2t — 50000(100 + 24)*/2. 
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Precisamos mais 300 litros de água salgada para o tanque ficar cheio. Mas a cada minuto entra 
2 litros de água salgada no tanque. Portanto o tanque estará cheio em 150 minutos. 


Isto é, quando t = 150. 


y = 100 + 2(150) — 50000(100 + 300) 3/2 
= 100 + 300 — 50000(400) 3/2 


= 400 — 50000(20)* 


= 400 - 6,25 = 393,75kg 


Portanto o tanque conterá 393,75 kg de sal quando estiver cheio de água salgada. 
1.7.2.1 Exercícios 


1. Um tanque contém inicialmente 100 litros de água fresca. Despeja-se então água contendo 
! E : : d ) 
3 kg de sal por litro no tanque, em uma taxa de 2 litros por minuto, e deixa-se a mistura 
escoar na mesma taxa. Após 10 minutos, interrompe-se o processo e adiciona-se água 
fresca numa taxa de 2 litros por minuto, com escoamento da mistura à mesma taxa. 


Determine a quantidade de sal no tanque ao fim de 20 minutos. 


2. Um tanque contém inicialmente 60 L de água pura. Água salgada contendo 1 kg de sal 
por litro entra no tanque a uma taxa de 2 L/min, e a solução (perfeitamente homogênea) 


deixa o tanque a uma taxa de 3 L/min; o tanque está vazio depois de 1 hora. 
(a) Ache a quantidade de sal no tanque depois de t minutos 


(b) Qual a quantidade máxima de sal alcançada no tanque? 


3. Um tanque contém inicialmente 100 L de água salgada com 1000 gramas de sal. Água 
salgada contendo uma concentração de sal de (1 + 0,05 sen t) gramas por litro entra no 
tanque a uma taxa de 2 L/min, e a solução (perfeitamente homogênea) deixa o tanque na 
mesma taxa. Seja y(t) a quantidade de sal no tanque no tempo t. 

Escreva uma equação diferencial para y e resolve ela para encontrar y explicitamente 


como uma função de t. 


4. Um tanque contém inicialmente 50 L de água salgada com 500 gramas de sal. Água 


salgada contendo 100 gramas de sal por litro entra no tanque a uma taxa de 2 L/min, e 
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a solução (perfeitamente homogênea) deixa o tanque a uma taxa de 6 L/min . Seja y(t) a 
quantidade de sal no tanque no tempo t. 
Escreva uma equação diferencial para y e resolve ela para encontrar y explicitamente 


como uma função de t. 


5. Um tanque contém 1000 L de água salgada com 15 kg de sal dissolvido. Água pura entra 
no tanque a uma taxa de 10 L/min. A solução é mantida bem misturada e sai do tanque 


na mesma taxa. Quanto sal permanece no tanque 
(a) depois de t minutos? 


(b) depois de 20 minutos? 


1.7.3 Modelagem de população 
Vamos modelar o fenômeno de crescimento da populações. 
Modelo simples 
e Seja P(t) o número de indivíduos em uma população no tempo t. 


e A população irá mudar com o tempo. De fato a taxa de variação de P será devido a 


natalidade (que aumenta P) e a mortalidade (que diminui P). Isto é 


a taxa de variação a taxa de a taxa de 


P no tempo t natalidade mortalidade 


e Vamos supor que todos os indivíduos são idênticos na população e que a taxa de nata- 
lidade média per capita, r, e a taxa de mortalidade média per capita, m, são constantes 


positivos. Isto é 


número de nascidos vivos por ano 


a 


—: ki 1 = 
r = taxa de natalidade per capita nanda populaão 


É É número de óbitos por ano 
m = taxa de mortalidade per capita = ————————— 
tamanho da população 


e Então o número total de nascidos vivos na população no ano t é rP, e o número total de 


mortes da população no ano t é mP. 
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dP 
e Portanto a taxa de variação da população, ET é dado por 


C=1P-mP=(r-mP=kp, 


onde k = (r — m). 


Portanto, a população cresce a uma taxa proporcional ao seu tamanho, isto satisfaz a equação 


diferencial 
dP 


ar = tP (1:72) 


onde k é chamado coeficiente de taxa de crescimento relativo, isto é, a taxa de natalidade menos 
a taxa de mortalidade. 


A equação ( 1.7.2) é de variável separável com solução 
P(t) = Po et. (1.7.3) 


Podemos deduzir as seguintes conclusões de (1.7.3): 


e Não estamos interessados em valores negativos de P, pois representa o tamanho da po- 


pulação. 


e a população irá crescer se k > 0 o que acontece quando r — m > 0 ou seja, quando a taxa 


de natalidade per capita, r, excede a taxa de mortalidade per capita 1. 


e Sek < 0, ou de forma equivalente, r < m então mais pessoas morrem em média do que 


nascem, de modo a que a população irá encolher e (eventualmente) ficar extinta. 


1.7.4 Problemas de crescimento e decaimento 


Seja a equação diferencial 


y =ky. (1.7.4) 


Se y = y(t), depende do tempo, esta equação diz que a taxa de variação de y é proporcional a y. 
Esta quantidade aumenta se k > 0 e diminui se k < 0. 


A variável y pode modelar: 


e o tamanho de uma população 
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e a quantidade de uma substância radioativa (neste caso k < O é a taxa de decaimento). 
e a quantidade de dinheiro investido. Neste caso, k é a taxa de juros. 
e etc. 


Como visto anteriormente, esta equação diferencial 1.7.4 é uma equação a variáveis sepa- 


ráveis cuja solução é y(t) = ce”. 


Exemplo 1.7.2. Calcule o rendimento do capital R$10.000,00 investidos por 5 anos a uma taxa 


de juros anual de 5%. 


Solução: A constante k = 0.05 e portanto, a solução é y = enter 


Como a condição inicial 
y(0) = 10000, segue que y(t) = 10000€ºº!. Após uma aplicação de 5 anos, o capital aumenta 
para: 


(5) = 10000€º02> = 12.840,25, 


isto é, o capital rendeu R$2.840,25. 


Exemplo 1.7.3. O Sudário de Turin é conhecido por ser o manto funerário de Jesus Cristo. Está 
preservado para exposição na Catedral San Giovanni Battista em Turin desde 1578. Três testes 
independentes de datação feitos em 1988 reveleram que a quantidade de carbono 14 do manto 
estaria entre 99,119% e 99,275% da quantidade encontrada em um tecido novo equivalente. 


Como podemos usar esses dados para determinar a idade do sudário? 


Solução: Este é um típico problema de decaimento radioativo, onde a taxa de decaimento (y”) 
é proporcional à quantidade do seu material (ky) a cada instante t. Para o carbono 14, sabe-se 
que a constante de proporcionalidade é k = —0,00001216. Assim, a equação diferencial que 
modela o problema é ' = —0,00001216y, cuja solução é y = ce “99901216! Chamamos yo = y(0) 
a quantidade inicial de carbono 14 no sudário e substituindo na solução, achamos o valor da 


constante c: 


—0,00001216x0 


y(O) = ce =cC=4y0 


isto é, a solução é y = yoe “2Zt6' Como dito no enunciado, a quantidade de carbono 14 em 
1988 era de 99,119% do valor inicial yo. Quanto tempo levou para a quantidade de carbono 
14 decair de yo para 0.99119y9? Substituindo estes dados na equação e resolvendo-a para a 
variável tempo t: 


0.99119yo = A > t=727,71anos, 
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isto é, o sudário deve ter sido fabricado por volta do ano 1260. Desta forma, fica provado que 


nunca poderia ser usado como manto funerário de Jesus Cristo. 
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UNIDADE 2 


EQUAÇÕES DIFERENCIAIS DE 
ORDEM SUPERIORES E A 
TRANSFORMADA DE LAPLACE 


Aula 2.1 


Equações Diferenciais de segunda 


ordem 





Exemplo 2.1.1. A equação diferencial y" = 2xyy' + e* é uma equação de segunda ordem, 


porém não é linear por causa do termo 2xyy”. 


Exemplo 2.1.2. A equação diferencial y” = 2xy' + e*y é uma equação de segunda ordem 
Pp q y bo, Ejtral 8 


linear. 
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Exemplo 2.1.3. A equação diferencial y” = 2xy + e * é uma equação de segunda ordem linear. 





Exemplo 2.1.4. A função (x) = e é solução da equação y" — 1/2y — 1/2y = 0. De fato, a 


função (x) é infinitamente diferenciável e ela satisfaz a equação pois e — 1/2 — 1/2 =0. 


Exemplo 2.1.5. A função q(x) = x definida em (0, +00) é solução da equação 
1 5 
ZA n2 
us +5>y. 
VE Ga) toy 


De fato, a função (x) é infinitamente diferenciável. Vamos verificar se ela satisfaz a equação 


diferencial: 


Í 
6x = (3x7)? + o > 6x=x+5x > í6x=6%. 
9x3 x2 





Exemplo 2.1.6. Seja o seguinte PVI 


W" +5y +6x=0 
y(O) = 2 
y(0) =3. 
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—2X 7 p3% 


Verifique que a função q(x) = 9e é solução deste PVI. De fato, a função q é pelo menos 


duas vezes diferenciável e satisfaz a equação diferencial: 


360 — 630% 900 2 +1050 +54 2 490% =) > 0=0. 


Além disso, a função q satisfaz as condições iniciais: (0) =9-7=2eqW'(0) =-18+21=3. 


Estudaremos métodos para resolver a equação (2.1.1) somente quando p(x) e g(x) forem 
constantes. Estes métodos poderão ser generalizados para equações diferenciais lineares de 
ordem superior. 

Como não é possível achar uma solução da equação (2.1.1) de forma simples e/ou apenas 
por inspeção, vamos estudar as condições que garantem que a equação (2.1.1) tenha uma 


solução. 


Teorema 2.1.1. [Teorema de Existência e Unicidade] Seja o problema de valor inicial: 


y" + pldy + ql)y = (x). 
Y(xo) = yo 


y (xo) = y1. 


Considere p(x), q(x) e r(x) funções contínuas num intervalo [a,b]. Se xo E [a,b] e se yo e yW forem 
números reais quaisquer, então a equação (2.1.1) tem uma e somente uma solução y(x) no intervalo 


[a, b]. 





Exemplo 2.1.7. Dado o problema de valor inicial 


ga "+xy= l 
y a = 


1 1 
onde as funções p(x) = —, q(x) = x e r(x) = — estão definidas no domínio [1,5]. Então, como 
P x? x2 


as funções p(x), q(x) e r(x) são contínuas no seu domínio e 2 E [1,5], o teorema de existência e 


unicidade nos garante que existe UMA e SOMENTE UMA solução y(x) do problema de valor 
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inicial definida no intervalo [1,5]. 


Exemplo 2.1.8. Mostre que y = x? sen x e y = 0 são ambas soluções de 
dy" -4xy +(2+6)y=0, 


e que ambas satisfazem as condições y(0) = 0 e y'(0) = O. Será que isto contradiz o Teorema 
da Existência e Unicidade de soluções? Note que as condições iniciais são calculadas no ponto 


x = 0 eneste ponto as funções p(x), q(x) er(x) não estão definidas. Logo, o teorema não se aplica. 


2.1.1 Exercícios 
1. Determine o intervalo de maior amplitude dentro do qual o problema de valor inicial 
proposto tem uma única solução. Não encontre a solução. 
(a) xy” +3y =x, y(l) = 1 e y'(1) = 0. Resp: JO, +oo[. 
(b) (x—- 1D)y” -3xy +4y= senx, y(-2)=2e y(-2) = 1. Resp: ] — 00,1[. 
(c) x(x — By” +5xy' +4y=3, y(3) = 0e y'(3) = —1. Resp: ]0, 4]. 


(d) (x—-3)y” + xy + (Inlxb)y = 0, y(1) =1 e y(1) = 0. Resp: ]0, 3[. 
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Aula 2.2 


Equações diferenciais lineares 
homogêneas de segunda ordem - 


Solução geral 


Estamos interessados em encontrar a solução geral da equação diferencial linear homogênea 
de segunda ordem: 


W +ploy +q0)y=0, (2.2.1) 
Sabemos da Aula 2.1 que um problema de valor inicial com a equação (2.2.1) contém duas 
condições iniciais: 


W+pody +aG)y=0, uxo)=V, víxo)=mn 


Neste caso, uma solução geral deve envolver duas constantes arbitrários. De fato, este é o caso. 
Mais ainda, vamos ver que existem duas funções diferentes que podemos combinar para criar 
uma solução para uma equação linear de segunda ordem. A técnica dependerá do princípio de 


superposição 
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2.2.1 Princípio da Superposição 


Teorema 2.2.1. [Princípio da Superposição] Se yi(x) e y>(x) são duas soluções da equação 


homogênea (2.2.1), então 


yO) = cry) + coyo(x) (2.2.2) 


é também uma solução para quaisquer constantes reais cy e co. 





Demonstração. Como y(x) e y>(x) são soluções de ( 2.2.1) temos que 


w +pl)y +q6)n =0 (1) 
Yo +pl)y +q(x)yo =0 (2) 


Substituímos y = ciyi(x) + c>y>(x) na equação homogênea: 


(cinto + cavalo) + plo(enn (o + cano) + q0Nenyn() + oyo) = 


cy + coy> + pl(ciy + coy5) + q((ciym + coy2) 


ci (vi + pldy + gy) +c> (y5 + pl)ys + qg(x)y>) 
Pq A iteração a o 
-0 de (1) =0 de (2) 
cqx0+cx0 


= 0. 


Portanto y = ciyi(x) + c>y>(x) é uma solução de (2.2.1). E 


Da Álgebra Linear, a solução ( 2.2.2) é chamada de combinação linear das soluções y1(x) e 
y>(x). O princípio da superposição diz que qualquer combinação linear de duas soluções da 
equação homogênea é também solução. 


Observação. Dadas duas soluções y e y> da equação 


y" + pld)y + qld)y =0, (221) 
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podemos construir infinitas outras soluções da forma 


y(x) = ciyi(x) + coyolx) (2.2.2) 


fazendo c4 e co percorrer todo o conjunto RR, isto é, 


y(x) = cryi(x) + coyolx) 


tem a forma da solução geral de (2.2.1). 


Portanto, uma pergunta surge: “Se y e y> são soluções de (2.2.1), será que a expressão 
(2.2.2) é a solução geral? Isto é, será que todas as soluções de ( 2.2.1) podem ser escritos como 
uma combinação linear de y e y> ou algumas soluções têm uma forma totalmente diferente”? 
Veremos a seguir que a equação (2.2.2) será solução geral se as soluções y e > satisfizerem 


uma determinada propriedade. 


Definição 2.2.1. Dizemos que duas soluções y e y> formam um conjunto fundamental de 
soluções de (2.2.1) se qualquer solução de (2.2.1) pode ser expressa como uma combinação 


linear de y e y>. 


Teorema 2.2.2 (Teorema da solução geral). Se p e q são funções continuas no intervalo aberto 


I=(a,b)esey e y> são duas soluções da equação diferencial linear homogênea: 


yW” + pl)y +gl)y =0 


satisfazendo 


Wíyi, yo)(xo) = yi(xo)yo(xo) — yi(xo)yalxo) + O 


para algum ponto xo € 1, então qualquer outra solução de (2.2.1) no intervalo I pode ser escrito 


unicamente da forma 


v=cym + Coy. 


A expressão y = cry + c>y> é chamado de solução geral da equação diferencial linear de segunda 


ordem homogênea. 
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Demonstração. Sejam y(x) e y>(x) duas soluções da equação (2.2.1) no intervalo I e seja Y uma 
outra solução em 1. Escolha um ponto xo € 1. Pela Teorema de Existência e Unicidade (Teorema 


2.1.1) sabemos que existe uma e somente uma solução y(x) da equação ( 2.2.1) tal que 


y(xo) = Y(xo) 
y (xo) = Y'(x0) 


(258) 


Portanto se pudermos mostrar que a solução da forma 


y(x) = crya(x) + coyo(x) 


satisfaz as condições iniciais (2.2.3), então teremos que Y(x) = ciy(x)+c>y>(x) é um combinação 
linear de y;(x) e y>(x). Portanto queremos mostrar que é possível escolher constantes c; e c> tais 
que: 


ciyilxo) + coyalxo) = (xo) =YVo 


y (xo) = Yo 


ciyi(xo) + coyo(xo) 


Para encontrar as constantes cj e c», precisamos resolver o sistema. Usaremos o método de 


Cramer: 


vo Y2(X0) vi(xo) vo 
vo Yo(xo) vi(xo) Yo 
g=————————— e g=——— 
vi(xo) ya(xo) yi(xo) Yya(xo) 
vi(xo) Yo(xo) vilxo) Yolxo) 


Portanto podemos resolver por cy e c> sempre que o denominador 


ias O e a a O 
vi(xo) Yolxo) 


Portanto sempre podemos expressa qualquer solução de (2.2.1) como uma combinação linear 


dey ey> se W(y,y>) £ 0. E 
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Exemplo 2.2.1. Considere a equação diferencial 


2%y' +3xy -y=0, x>0 


1/2 


Mostre que as funções y = x!2 e y =x! formam um conjunto fundamental de soluções e 


escreva a solução geral. 
Solução: 


e ye y> são soluções pois 








— "3/2 —1/2 
Ze 43x yr = 22 [ o Jea(5 J-ue=(Se5-a)tt-o 


222y) + 3xy) — yo = 2x2 (2º) + 3x(-x?) = =(4-3-=1)x1=0 


e Para mostrar que y e y> formam um conjunto fundamental de soluções, calculamos o 
Wronskiano de y e yo : 
1/2 -1 


yWw 4Yy2 x 


A Y% 


1 3 -3 
= = = -y3/2. 2y 32 Ly 
. 2" o 2448 














1 ap 
2* X 


Como W £ O parax>0, y, y> forma um conjunto fundamental de soluções e a solução geral 
é dada por 


C2 
yY=ciy +Coyo = cx? + Fi 


onde c4, c são constantes. 
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2.2.2 Wronskiano e Independência Linear 


Definição 2.2.3 (Independência linear). Dizemos que duas funções reais yi(x) e y>(x) são line- 
armente independentes no intervalo (a, b) se não é possível escrever uma delas como múltiplo 
da outra, isto é, não existe c E R que satisfaça a relação y(x) = cy>(x). Se duas funções não são 


linearmente independentes, dizemos que elas são linearmente dependentes. 


Observação. Outra maneira de demostrar dependência linear é provar que podemos ter 


ciyi(x) + coyolx) = O 


para todo x E (a,b) com cy e c> não ambos nulos. 


2 
x, x>0 

Exemplo 2.2.2. Seja y(x) = xº e y>(x) = xx = E: - Mostre que y e y> são 
—xº, x<0 


linearmente independentes em R. 
Solução: Suponha que y; e y> são linearmente dependentes em RR, isto é 
2 
ciyi(x) + coyo(x) = cyx” + cox)x] = 0 
sendo que c; e c, não são ambos zeros. Colocando x = 1 e x = —1 temos o sistema 


c+0c=0 


c—c=0 


Resolvendo o sistema temos que c; = c> = 0, uma contradição à nossa hipótese. Portanto yj e 
y> são linearmente independentes em R. 


O próximo resultado carateriza independência linear em termos do Wronskiano. 
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Teorema 2.2.3. Sejam y e y> funções deriváveis no intervalo aberto (a, b). 


(1) Se y e y> são linearmente dependentes em (a, b) então W(y, yo)(x) =0 Yx e (a,b) 


(2) Se W(y, y>)(xo) £ O para algum xo E (a,b) então y e y> são linearmente independentes em 


(a, b). 





Demonstração. Observe que a segunda afirmação é a contra-positiva da primeira, então vamos 
demonstrar a primeira. 

Suponha que as funções y e y são linearmente dependentes no (a, b). Então para todo 
x € (a, b) temos que 


yi(x) = cyo(x) 


Portanto o Wronskiano é 


cya(x) yo(x) 
cyslx) yo(x) 


yi(x) yo(x) 
vil) yo(x) 


Wíy, Yy2) = = cya(x)ys(x) — cys(x)yo(x) = 0. 














O próximo resultado é muito importante pois ele dá uma fórmula simples para o Wronski- 


ano. 


Teorema 2.2.4 (Fórmula de Abel). Sejam p e q funções contínuas em um intervalo (a, b). Sejam 


y1 e y> soluções da equação diferencial 


y" + pldy” + gl)y = 0. 


Então o Wronskiano W(w1, y>)(x) é dado pela fórmula 


WC, yo)(x) =: e Srt) dx 
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Demonstração. Pela definição, 
W=nymp-AnrS W=yy+nym-ywmy-nhymS W=inyy-y( yo 
Como y e y> são soluções de (2.2.1), então: 


vi (00) +pld)y + gl) =0 eva (x) + pl)ys + glx)yo = 0. 


" 


vo) = py - qn ey(x) =-plo)y, — q(x)yo. 


tn 


Substituindo y/ e y; na expressão de W”, obtemos: 


, 


Wº = n(-pl)yo — q6)yo) = yal-pody — q6)yn). 
Colocando p(x) e q(x) em evidência na equação acima, temos: 
Wº = pl yo + voy) + 0N-yny +inyo) > W=-pW 
Resolvendo a equação diferencial para o wronskiano W, temos 


W=-9W > W(6)= cem pedax, 


Exemplo 2.2.3. Calcule o Wronskiano das duas soluções da equação 


sen E) 


dy" -xy +(cos?x+e v=0, x>0. 


Solução: Como x > 0, podemos reescrever a equação diferencial na forma padrão 
7, 1 , 1 2 sen x) — 
V=y += (cos XE )=0 
1 j 
Portanto p(x) = E Logo pela Fórmula de Abel temos que 


Nise lradr= des par 
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Teorema 2.2.5 (Teorema do Wronkiano). Se yi(x) e y>(x) são duas soluções da equação 


homogênea (2.2.1) definidas em (a, b), então seu Wronskiano é: 


Wíy, y2)(x) =0, Vx € (a, b). 


Wíy,yo)lx) £O, Vxe (a,b). 





Demonstração. Pela Fórmula de Abel temos que 
Wí(x) = c- e) r(o da 


e Sec=0,entãoW=0,Yxe (a,b). 


e Caso contrário, sec + 0, W £ 0, Yx E (a,b), pois W é uma exponencial e esta nunca se 


anula. 
E 


Para resumir a teoria, podemos relacionar todos os resultados sobre conjunto fundamental 
de soluções, wronskianos e independência linear da seguinte forma: sejam y e y> soluções 
da equação homogênea 


yr + pe)y + q6)y =0. 
As seguintes quatro afirmações são equivalente: 
e As funções 1 e y> formam um conjunto fundamental de soluções no intervalo T. 
e As funções 11 e y> são linearmente independentes em T. 
e W(y, y>)(xo) £ 0, para algum xg € 1. 
e Wí(y, y>)(x) £ 0, para todo x e T. 


2.2.3 Exercícios 


1. Calcule o Wronskiano dos seguintes pares de funções: 
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(a) er qe Resp: Se! EA 

(b) cost, sent. Resp: 1. 

(c) Rd Resp: et. 

(d) e sent, e cost. Resp: e? 


(e) cos 28, 1+ cos 20. Resp: 0. 


. Verifique que y1(t) = É e y>(t) = t! são duas soluções da equação diferencial É y' —-2y = 0 


para t > 0. Então mostre que y(t) = cit? + c»t! é também uma solução desta equação para 


quaisquer cj e c». 


- Mostre que y(x) = cy sen x+c> cos x é uma solução geral de y' +y = Oem algumintervalo. 


A seguir, ache a solução do PVI para o qual y(0) = 2e y'(0) = 3. Resp: y = 3 sen x+2 cos x. 


?2 são soluções de yy” + (y)? = 0, para x > 0. Mostrar 


então que cy + cox!/? 


não é em geral solução desta equação. Por que não? Isso contradiz o 
teorema da existência e unicidade de soluções para equações lineares de segunda ordem? 


Resp: A equação não é linear. 


. Seo wronskiano W de fe g é 3” e f(x) = e”, ache g(x). Resp: 3xe”* + ce? 
. Seo wronskiano W de fe g é x2e" e f(x) = x, ache g(x). Resp: xe! + cx. 


- Se W(f,8) é o wronskiano de fegeseu=2f -gev= f+2g9,ache o wronskiano W(u,v) 


de ue v em termos do wronkiano W(f, 2) de fe g. Resp: 5W(f, 9). 


- SeW(f,g)=xcosx- senx, eseu=f+3gev=f-g,acheo wronskiano W(u,v) deu e 


v. Resp: —4(x cos x — sen 2). 


- Em cada uma das alíneas abaixo, verifique que as funções dadas y e y> são soluções da 


equação diferencial. Elas formam um conjunto fundamental de soluções? 


(a) y" +4y=0,yi(x) = cos 2x, y(x) = sen 2x. Resp: Sim. 
(b) y"-2y +y=0, 1n(x) = e, yolx) = xe”. Resp: Sim. 
(9) 2y" xx +Dy +(x+2)y=0, x>0,1n(x) =x, y>(x) = xe”. Resp: Sim. 


O wronkiano de duas funções é W(x) = x sen ? 


x. As duas funções são linearmente 
dependentes ou independentes? Justifique. 


Resp: São linearmente independentes porque o wronkiano não se anula em todo ponto. 
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11. 


12. 


13. 


14. 


Ache o wronskiano das duas soluções da equação diferencial dada sem resolvê-la. 


(a) x2y” — x(x + 2)y + (x + 2)y = 0. Resp: cx?e”. 


(b) x2y” + xy' + (12 — 4)y = 0. Resp: c/x 


Mostre que se p(x) é diferenciável e p(x) > 0, então o wronskiano de duas soluções de 


[py] + q(x)y = O é W(x) = c/p(x), onde c é uma constante. 


Se y e y são duas soluções linearmente independentes de xy” + 2y' + xe'y = 0 e se 


W(y, y>)(1) = 2, ache o valor de W(yr, y>)(5). Resp: 2/25 


Se y1 e y> são duas soluções linearmente independentes de 22" —2y +(3+x)y=0ese 


W(1n, y>)(2) = 3, ache o valor de W(y,, y>)(4). Resp: 3 ve. 
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Aula 2.3 


O Método de Redução de Ordem 


O Método de Redução de ordem é um método para converter uma equação diferencial 
linear para uma equação diferencial linear de ordem inferior, e, em seguida construir a solução 
geral da equação diferencial original usando a solução geral da equação de inferior. 

Vamos ilustrar o método de "redução de ordem"para equações diferenciais lineares homo- 
gêneas de segunda ordem. 


Problema: Seja a equação diferencial linear homogênea de ordem 2. 


yw” + pl)y + q(u)y =0, (2.3.1) 


onde p(x) e g(x) são quaisquer funções contínuas definidas no intervalo (a,b). Dadas duas 


soluções y e y> linearmente independentes desta equação, é fácil escrever a sua solução geral: 


y(x) = cryi(x) + yo(x). 


Porém, infelizmente não existe um método geral para encontrar estas soluções LI da equação 
diferencial (2.3.1) no seu caso geral. 

Suponha, no entanto, que conhecemos uma das soluções y da equação homogênea. O 
que podemos fazer para encontrar uma segunda solução > tal que o conjunto [y1, y>) seja um 
conjunto fundamental de soluções? 


Apresentamos dois métodos: 


2.3.1 Método 1: Usando a Fórmula de Abel 


Neste método escrevemos o Wronskiano W[11, >] em 2 maneiras: 
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(1) Pela definição: WlIy, y2] = yiy5 — y2W 
(2) Pela Fórmula de Abel: Wly,, 12] = ce” Seo dx 


Igualando as duas teremos uma equação diferencial de primeira ordem de y> que podemos 


resolver. Veremos este procedimento no seguinte exemplo: 


Exemplo 2.3.1. Seja y = x uma solução da equação diferencial homogênea de segunda ordem 
2,4 


Xu +xy -4=0, 2>0. 


Encontre a a sua solução geral. 
Solução: Primeiramente, dividimos a equação diferencial por x?: 


11 l , ú E 
eae 


Seja y> a outra solução linearmente independente com 11. 


e Pela definição do Wronskiano 
Wiy — 1, y2] = y1y> — Ya4h = xy — 42 
e Pela Fórmula de Abel 
Wlyi, y2] = ce JPG) dx — co J(G) dr — corlnx — - = - (fazendo c = 1) 
Igualando os dois valores do Wronskiano temos a equação diferencial linear de 1º ordem 


: l | 1 
pp hocçhe= a 


Equação é linear. 
1 
Fator integrante: e Sad = erlnx = 


Multiplicando a equação pelo fator integrante temos 
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Portanto 


Integrando temos 


pm 
te 
[9] 

Il 
q 
L 
ou 
R 


1 1 

e = o 
1 

un =x 


Portanto, a solução geral é: 


1 
= = + — 
RAS 2x 


2.3.2 Método 2: D'Alembert 


Seja y1(x) uma solução de (2.3.1) diferente de zero e assim, cy; é também uma solução de 


(2.3.1) para qualquer constante c:. 

Procuramos a segunda solução y> de tal forma a garantir a independência linear das duas 
soluções y1 e 12, isto é, que y>/y1 + constante. Daí, existe uma função v(x) tal que y>/y1 = v(x) e 
a solução 1» será escrita na forma: 


Ya() = (xo) ya (29). 


Precisamos então determinar v(x) de tal forma que y> = vy; seja solução da equação (2.3.1). 


Veremos este procedimento no seguinte exemplo: 


Exemplo 2.3.2. Seja y = x uma solução da equação diferencial homogênea de segunda ordem 


Cy +xy -y=0, x>0 


Encontre a a sua solução geral. 


Solução: Primeiramente, dividimos a equação diferencial por x2: y” + (1/x)y — (1/x2)y = 0. 
Fazendo y> = vy;, uma solução da equação diferencial, substituímos > e suas derivadas na 
equação: 


p=zn5yp=vn+vyy=pn+rovy+oy +vy,. 
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Substituindo y = x,y, = 1 ey] = O na equação acima, temos: 
p=uSyW=vrto>sy =v'x+20'x. 


Finalmente, substituindo agora y> e suas derivadas na equação: diferencial y” + (1/x)y' — 


(1/22)y = 0, temos: 


1 1 
vx+2 + (vx+0)- Sux=0 = x" 420 40 420-2=0 
é x ua 
TAVA TAVA 3 
= ia do = RE aa = [Sar [Sa 5 
x v x v x 
> Inv =Inx”? Em => as A > ai ad 
. o 22 = 


1 
Logo, a solução geral é y = cix + 5 


2.3.3 Exercícios 


1. Em cada problema abaixo, ache uma segunda solução da equação diferencial dada, co- 


nhecendo uma solução. 


(a) 2y" +3xy -y=0,x>0,  ()=x!. Resp: = ma 
(b) («-Dy'-xy/ +y=0,x>1,  (g)=e. Resp:y=x. 


(ow +5y=0 mwm=l. Resp: y = se 


(dy -/=0 mn=1. Respy=e 
(e) y" — 4y' +4y=0; n=e. Resp: y = xe” 


() y' +16y=0; wi = cos 4x. Resp: y= sen 4x 
In(1 + e9) + ! 


(8) y' —-y=0; y1 = coshx. Resp: y = coshx( TE 


(h) 9y” — 12y' +4y=0; =" Resp: y = xe3* 
4 


(1) 2y” - 7xy' + 16y=0; Wn=%. Resp: y=x*Inx 


pp (1-2x-y" +21 +0)y -24=0; yw=x+1. Resp: y=2)+x+2 
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Aula 2.4 


Equações diferenciais lineares 
homogêneas de coeficientes constantes 


de segunda ordem 


Vamos agora estudar as equações homogêneas ( 2.3.1) com coeficientes constantes, isto é, 


p(x) = a e q(x) = b, ondea e b são constantes: 
y"+ay +by=0. (2.4.1) 


Observe que a função exponencial e** tem a propriedade de que suas derivadas são todas 


A 


constantes multiplicadas por si própria. Isto nos faz considerar a função y = e** como uma 


possível solução de ( 2.4.1) se À por convenientemente escolhida. 
2.4.1 Determinando o valor de À 


Fazendo y = e” solução da equação e substituindo y e suas derivadas y' = Ae" ey” = AZe** 


na equação diferencial (2.4.1), temos: 


y+ay +by = 0 
= Ae ra shot = 0 
se +ar+b) = 0 

> Mrad+b = 0. (2.4.2) 
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Observe que a função exponencial e” nunca se anula. A equação (2.4.2) é chamada de 
equação característica associada à equação diferencial (2.4.1). Para acharmos o valor de 4, 


basta resolver a equação característica. 


a E am <= = 2 — 
Ee a pas (2.4.3) 


Anis 2 2 


Podemos ter três casos para os valores de 4: 
e Caso 1: as raízes são reais e distintas. 
e Caso 2: as raízes são complexas. 
e Caso 3: as raízes são reais e iguais 


2.4.1.1 Caso 1: raízes reais e distintas 


Se as raízes da equação característica A, e À» são reais e distintas, as soluções da equação 
diferencial serão: 


n=d” e p=" (2.4.4) 


Observação: É fácil ver que estas duas soluções y1 e y> são linearmente independentes uma 


vez que 


Aox 
E = — = eo A) * constante. 
1 € 
Logo, a solução geral é 
u(x) = cre + coefe (2.4.5) 


Exemplo 2.4.1. Resolva o seguinte problema de valor inicial: 


W+y-2y=0, 


v(O) =1 
y(0)=1 
Solução: A equação característica é 2 + À — 2 = Q cuja solução é: 4 = 1e Ap = -2. Desta 


forma, a solução geral é: 


y(x) = ce" + coe 
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Substituindo as condições iniciais, temos: 


q+0o=1, 


c1—20)=1 


Resolvendo o sistema acima, obtemos: c; = 1e c) = 0 e a solução do PVI é: y(x) = e”. 


2.4.1.2 Caso 2: raízes complexas conjugadas 


Neste caso, as raízes podem ser escritas na forma: 
M=r+10 e Ao =r-—10, 
onde 1 é o número imaginário. Assim, 


elrtio)x — 


e% pix = e ( 


cos 0x + i sen 0x). 


Note que a expressão e? = cos 0 +i sen O é conhecida como fórmula de Euler. De forma 


análoga: 


elr-iO)x = ex gTiOx = elx ( cos Ox -— 1 sen ga). 


Como estamos interessados apenas em soluções reais e 


elHiOx = EX ( cos Ox + i sen 0x) 


eli = EX ( cos Ox — i sen 67), 


então podemos obter 1 e y> da seguinte forma: 


(r+i0)x (r-i0)x 
e +e 
n="20 E e” cos 0x. 


gli) Ee elr-io)x 
RE. e” sen 0x. 


e a solução geral é: 


y(x) = ce” cos Ox + coe” sen Ox. (2.4.6) 


Joseph N. A. Yartey & Simone S. Ribeiro 


98 


Exemplo 2.4.2. Resolva o seguinte problema de valor inicial: 


2y" -4y' +8y=0, 
y(O) = 1 
y(0) = 0 


Solução: A equação característica é A2-2A4+4=0 cuja solução é: 4 = 1 + V3ieA, =1- V3i 


Neste caso, a solução geral é: 
y(x) = cj" cos V3x + coº sen V3x. 
Substituindo as condições iniciais, temos: 


cg=1, 


q +c2 V3 =1 
Resolvendo o sistema acima, obtemos: c = 1ec) =-—1/ vB. Daí, a solução do PVI é: 


1 
(x) = e" cos V3x—- —e sen V3x. 
l 5 


2.4.1.3 Caso 3: raízes reais e iguais 


Neste caso, a única raiz repetida da equação característica é À = —a/2 e assim, teremos 
apenas uma solução: 


Vs ea, 


Precisamos determinar uma segunda solução y> que seja linearmente independente com a 


primeira. Como já vimos antes, fazemos y> = v(x)y; e encontramos a segunda solução. 


Exercício 2.4.1. Prove que a segunda solução é yp = xe "*!2, 


Exemplo 2.4.3. Resolva o seguinte problema de valor inicial: 


y" +4y +4y=0, 
Y(O) = —2 
y(0) =7 
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Solução: A equação característica é A? + 4) + 4 = 0 cuja solução é: À = —2 com multiplicidade 
2. Daí, a solução geral é: 


y(x) = ce + coxo. 


Substituindo as condições iniciais, temos: 


Cc — —2, 


—2c +c0=7 


Resolvendo o sistema acima, obtemos: cy = -2 ec, = 3. Daí, a solução do PVI é: 


u(x) = 20 2 + 3xe 2 


2.4.2 Equação de Euler-Cauchy de segunda ordem 


As equações de Euler-Cauchy de segunda ordem são equações diferenciais do tipo: 


o2y" +mxy +agy=0 (2.4.7) 


No ponto x = 0, o termo a,x?y” da equação (2.4.7) se anula e por esse razão dizemos que 
x = 0 é um ponto singular para a equação de Euler-Cauchy. Qualquer solução de ( 2.4.7) estará 


definida para x > 0 ou para x < 0. 


O método para resolver a equação de Euler-Cauchy tem como principio a mudança de 


variável 
x=e 


que transforma a equação numa equação diferencial linear com coeficientes constantes nas 


variáveis y e t. 


e Pela regra da cadeia, 
E 
dr” dx dr “ax 
ou 
dy 1dy 


ERR (2.4.8) 
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e 
dy  dafdy) d Ay 
a? — arlar)” ala 
RR Po 
“dxl dx) di 
dy dy py dy 
E [+ E) = 2º qu "ax 
ou ; ; 
dy 1 /dy dy 
Substituindo equações (11.8) e (11.9) em (2.4.7) temos 
dy dy E 
va + (ay — 2) + a0y = 0. (2.4.10) 


que é uma equação diferencial de coeficientes constantes. 


Exemplo 2.4.4. Resolva 


d? d 
a no +2y=0. 


Solução: A equação é do tipo Euler-Cauchy, logo fazemos a mudança x = é out = Inx,e 


transformará a equação para 
dy 3 dy 


q? a +2y=0 


Egn. característica: 


A -3A+2=0>A=1, A2=2 
Portanto a solução geral é 


v=ce +cel =cx+cx 


pois x = é. 


24.3 Exercícios 
1. Ache a solução geral de cada uma das equações diferenciais abaixo: 


(a) y'+2y -3y=0. Resp: y=ce + ce. 
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(b) y'-2y +10y=0. Resp: y=cie' cos 3t + coe' sen 3t. 


(0) 64" -y -y=0. Resp: y=ce!? + ce 3. 


3 3 


(d) y" +6y' +13y=0. Resp: y=c e“ cos 2x+ ce" sen 2x. 


X 


(e) y'+5y =0. Resp: y=ci + ce. 


/2 


(f) 24" +2y +y=0. Resp: y=c e”? cos (x/2) + coe *!2 sen (x/2). 


(g) 4y" +9y =0. Resp: y=c; cos (3x/2) + co sen (3x/2). 
2. Ache a solução dos seguintes problemas de valor inicial: 
(a) y +y -2y=0,y(0)=1, y'(0) = 1. 
Resp: y=€ 
(b) y” —-6y +9y=0,y(0) = 0, y'(0) = 2. 
Resp: y = 2xe*. 
()yW+y=0, ,y(n/3)=2ey'(m/3) = —4. 
Resp: (1 +2 43) cos x— (2— 3) sen x. 
(d) y" +3y' =0, y(0) = —2, y'(0) = 3. 
Resp: y=-1- e, 


(e) y' +8y -9y=0,y(1)=1,y(1) =0. 


1 9 
R «yv=— —9(x—-1) 4 À a 
esp: Y= 9º + 19º 
() 9y” — 12y' + 4y =0, y(0) =2, yí(0) = —1. 
Resp: y = 2e*/8 — axei 
(9) y'+y+1.25y=0,y(0)=3, y'(0) = 1. 
Resp: y = 3€ “2 cos x + ae! 2 sen x 


3. Ache a equação diferencial cuja solução geral é y(x) = cj! + pe”?! 
Resp: y' +y —6y=0. 
4. Ache a equação diferencial cuja solução geral é y(x) = cje !!? + cpe 2. 


Resp: 2y" +5y' +2y=0. 


5. Considere o problema de valor inicial: 


4)" +12y' +9y=0, y0)=1, y(0)=-4. 
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> 
(a) Ache a solução deste problema. Resp: y = €2%/2 — sue a 


(b) Determine onde a solução é zero. Resp: x = 2/5. 


(c) Determine as coordenadas (xo, yo) do ponto de mínimo.Resp: xo = 16/15 e yo = 
Dor8I5 

3 à 

(d) Mude a segunda condição inicial para y'(0) = b e ache a solução como uma função 


de b. A seguir, ache o valor crítico de b que separa a solução na sua parte positiva da 


parte que pode ser negativa. Resp: y = e 2/2 + (b + 3/2)xe **/2 


6. Ache a solução do problema de valor inicial 


1 ai E > é = = 

pPeyso 0 Mb 

Calcule o seu valor mínimo no intervalo O < t <2. 

1 
Resp: y = zo +e*. Mínimo éy=1lemt=In2. 
7. Ache a solução do problema de valor inicial 
MA , " 1 
2y —-3y +y=0, 0 W0=2 y(0)==. 


Calcule o seu valor máximo e o ponto onde a função se anula. 


Resp: y = —e! + 3e!/2, Máximo é y = : emt=In(9/4)ey=0emt=In9. 
8. Considere o problema de valor inicial: 
4y" +4y' +y=0, y0)=1, y(0)=2. 


(a) Ache a solução deste problema. 


Resp: y = 2 + axe 


(b) Determine as coordenadas (to, yo) do ponto de máximo. 


Resp: xo = 8/5e yo = Se 415, 


(c) Mude a segunda condição inicial para y'(0) = b > 0 e ache a solução dependente de 
b. 


2b+1 
Resp: y=€*2 + xe *2, 
espiy=e 3 —*e 
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(d) Determine as coordenadas (tm, ym) do ponto de máximo em termos de b. Descreve 
como os valores de tm e ym quando b cresce. 


Resp: xm = ,ym = (1+2b)e PICA xy — 2 ym — 00 quando b — oo, 


4b 
2b+1 
9. Ache a solução do problema de valor inicial 


v -y =27=0, (O) = a, y'(0) = 2. 


Ache o valor de a para que a solução tenda a zero quando t tende a infinito. 


Resp: a = —2 
10. Ache a solução do problema de valor inicial 
4y" —- y=0, y(O) = 2, v'(0) = p. 


Ache o valor de 8 para que a solução tenda a zero quando t tende a infinito. 


Resp: B = —1. 


11. Considere o problema de valor inicial: 
yW +2y +6y=0, y0)=2, y(0)=a>o0. 


(a) Ache a solução deste problema. 


+2 
Resp: y=€ (2 cos V5x + E 





sen 52). 


(b) Ache a para que y = O quando t = 1. 
Resp: a = 1.5086. 
(c) Ache, como uma função de a, o menor valor positivo de x para o qual y = 0. 


245, 


qa+2 





1 
Resp: x = — arctg(— 
V5 
(d) Determine o limite da expressão da alínea (c) quando a — co. Resp: x = 7/ V5. 
12. Para cada uma das equações diferenciais abaixo, determine: 


(a) O valor de a, se existir, para o qual todas as soluções tendem a zero quando t —» oo. 


(b) O valor de a, se existir, para o qual todas as soluções não nulas se tornem ilimitadas 


quando t — oo. 
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iy”-(Qa-Dy +a(a — 1)y = 0. 
Resp: y — O para a < 0. yse torna ilimitado para a > 1. 


iy'+(3-oy -2a-1y=0. 


Resp: y — O para a < 1. Não existe a para o qual a solução se torna ilimitada. 


13. (Redução à primeira ordem.) Faça a substituição u = y' nas equações diferenciais abaixo 


e, a seguir, resolva-as. 


(a) x22y" +2xy —-1=0,x>0. 
Resp: y=cx!+c+lnx. 
(b) xy" +y' =1,x>0. 
Resp: y=cilnx+c +41. 
(O) y+uy7 =. 
Resp: y(x) = a arctg(x/c2) + c1. 
(D) y +y =e?. 
Resp: y(x) =e “(-c-x—1)+c, 


14. (a) Considere a equação y'+2ay +a?y = 0. Mostre que as raízes da equação característica 


são A = À2 = —a. Portanto, uma solução é y = e“. 


(b) A teoria do Wronskiano nos diz que 


Wíy, 42) = 14 — WHY = crer, 


onde cj é uma constante. 


(c) Use a alínea anterior para mostrar que a segunda solução é y>(x) = xe “*. 


15. Encontre a solução geral das seguintes equações de Euler-Cauchy: 


dy dy 
2 » RRDNEEDERS (7: 3/2 
(a) 4x q CEEE +3y=0 Resp: y=cx+cox 
dy dy 
2 Esq 7 
(b) x Tê — 6x am =0 Resp: y=C +c2x 
ade d 
(c) x A - 368 +4y=0 Resp: y= (cg +coInx)x? 
= Pp: y 
dx? dx 
pdy — dy 
— esp: y=cy cos (2ln 2 sen (2In 
(d) x “o tt +4y=0 Resp: y=c Qlnx)+c (21In x) 
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Aula 2.5 


Equações diferenciais lineares 


não-homogêneas de segunda ordem 


Considere a equação diferencial linear não-homogêneas de segunda ordem 
y" + pldy + q0)y = g() (2.5.1) 


onde g(x) £ 0. 


Teorema 2.5.1. Suponha que Yi(x) e Y>(x) são duas soluções de (2.5.1) então suas diferenças 
Y(2) = Ya (x) — Yo(x) 
é uma solução da equação homogênea correspondente 


y" + p(o)y" + a(x)y = 0. 


Se, além disso, y e y> forma um conjunto fundamental de soluções de ( 2.2.1), então 


Yi(x) — Yo(x) = cryilx) + coyo(x). 
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Demonstração. Vamos provar que Y(x) = Yi(x) — Y>(x) é uma solução de ( 2.2.1): 


Ye por + g(o)Y 


E emos =) cotas a 
(47 + pOr + g6)) = (15 + pe) + go) 


= 860)-g6)=0 


Portanto Y1(x) — Y>(x) é uma solução de ( 2.2.1) e como (1, y>) forma um conjunto fundamental 


de soluções de (2.2.1) temos que 


Ya (2x) — Yo(x) = cryalx) + coyo(x). 


Teorema 2.5.2 (Teorema da solução geral). Dado uma solução yp(x) da equação diferencial 
não-homogênea 


y" + pl)y' + ql)y = g(x) (2.5.1) 


então a solução geral da equação ( 2.5.1) é da forma 


yOo) = Yplx) + crya(x) + coyolx) (2.5.2) 


onde y1 e y> são duas soluções linearmente independentes da equação diferencial homogênea corres- 
pondente 


y" + pl)y + q(x)y = 0. (22.1) 





Demonstração. Sejam y(x) e yp(x) duas soluções de (2.5.1) então pelo Teorema 2.5.1 temos que 


y(x) — yplx) = cryi(x) + coyolx) 


Portanto 


y(X) = yp(x) + crya(x) + coyo(x) 
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Observação. Chamamos yc = ciyi(x) +c>y>(x) a solução complementar e yp(x) uma solução particular. 


Portanto a solução geral da equação (2.5.1) é da forma 
YO) = Vel) + Yp(x) 
Isto é, para determinar a solução geral da equação linear não-homogênea ( 2.5.1), seguimos os 3 passos: 
1. Determine a solução geral yc = ciyi(x) + coy>(x) da equação homogênea correspondente (2.2.1). 
2. Determine uma solução particular yp(x) da equação diferencial não-homogênea ( 2.5.1). 
3. Escreve a solução geral de ( 2.5.1) 


y(x) = ye) + yp(x) 


2.5.1 Soluções particulares 


As duas próximas seções vamos discutir duas métodos para determinar uma solução par- 


ticular y, da equação diferencial linear não homogênea de segunda ordem. 


2.5.1.1 Método de Variação de parâmetros 


Considere a equação diferencial linear não-homogêneas de segunda ordem 


yr + pO)y + q6)y = g(x) (2.5.1) 


onde g(x) £ 0. 
O método inicia-se supondo conhecidas duas soluções linearmente independentes y(x) e 


y>(x) da equação homogênea associada 
y” + pld)y + g(x)y = 0. (22.1) 
Hipótese: Procuramos uma solução particular da equação (2.5.1) da forma: 


Yo) = un (2) yr (x) + uol) yo (x) (2.5.3) 
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onde u1(x) e u>(x) são funções a determinar, satisfazendo a 


u6)yn(x) + us(a)yo(x) = O (2.5.4) 


Derivando (2.5.3) temos 


/ , , , 
Wy + Uiy + Usy2 + Uzy> 


a 
| 


uy + u2y, usando (2.5.4) 
Derivando y, temos 
A Pa MÁ Paul MÁ 
Yp = My +Uy + Usy + U2y 


Substituindo yp, Yp e Yp na equação (2.5.1) e efetuando algumas simplificações obtemos 


mnly) + pO)y + q0)ya] + uolys + pod + q()ya] + us + usyo = g(%) 


e como 1 e y> são soluções da ( 2.2.1), as expressões entre colchetes na última igualdade são 


iguais a zero e portanto 


Uy1 +Uoyo = 8(%) (2.5.5) 


Juntando as equações ( 2.5.4) e (2.5.5) temos o sistema linear nas variáveis uj e us : 


So 


my +usyo =0 mw y || 
ou 


uy + uy = g(%) vi dó || (3º) 


OQ 


Usando a Regra de Cramer deduzimos que 


; —y28(x) —ya g(%) 

grs Es a o (2.5.6) 
Comy-yiyo Wíiy,y) 

Sic = AB em ASR (2.57) 


2 ys-yiyo  W(y,yo) 


onde W(y:, y>) é o Wronskiano das funções 1 e y>, que por hipótese é diferente de zero. 


Podemos determinar as duas funções (x) e u2(x) integrando os resultados em (12.6) e (12.7). 
Resumo do Método: 


Para resolver y” + p()y' + q(x)y = g(x) 
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(1) Encontre o conjunto fundamental de soluções [y1, >] da equação homogênea associada. 


Portanto a solução complementar é yc = ciy1 + Coy> 
(2) Calcule o Wronskiano W(y, y>) 
(3) Determine 1 e u>, integrando os resultados em (12.6) e (12.7). 
(4) Uma solução particular é y, = uy + u2y>. 


Exemplo 2.5.1. Resolve a equação 
y"-3y +2y=e senx 
Solução: 


e Determine duas soluções 1.i na equação homogênea: 


Equação característica: A2-3A4+2=0> A =leA,=2. 


2. 


Portanto y = ee y> = e” e a função complementar é 


Ve= qe + ce 


e: e: 


a = 963% — 53X — 53x 
Ga a PR a 


e Observe que g(x) = e" sen x. Agora 


u = —yog(x) . -e* sen x dt 
t Wí(y, 42) ai | 


Integrando temos que m(x) = cos x. 


ns)  Xsenx E 


Dm 
Us = rr DN 


=X, 
WO, 95) Es e“ sena. 


—x 
Integrando temos que u2(x) = ea cos x + sen x). 
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e Uma solução particular é dada por 


Yp uy + Uay2 


ne 


e 
= e cosx- Eae cos x + sen 3)e” 


1 
= e cosx- 3º ( cos x + sen x) 


1 
= 3 ( cos x— sen x) 
Portanto a solução é 


1 
v=0ó+ cera 3º ( cos x— sen x) 


Exemplo 2.5.2. Use o método de variação de parâmetros para achar a solução geral da equação 
de Euler-Cauchy 


2y' +xy +y= sec(lnx), x>0. 


Solução: Equação é do tipo Euler-Cauchy, logo fazendo a mudança x = é, transformará a 
equação para a forma 

E () (2.5.8) 
— = sec Da 
de 
e Determine duas soluções L.I da equação homogênea associada a equação (12.8). 


Equação característica: 22 +1 = 0 > À = +i. Portanto as duas soluções L.I são 


yWw= cost e y= sent 


t Ê 
e W(y,y2) = Poa E (9) en = cos*(t)+ sen) =1 


vw yW — sen (t) cos (t) 


e Observe que g(t) = sec (t). Agora 


, o —ystt)  —sen(bsec(t) 
Fa WGy yo) 1 == teto 


Integrando temos que m(t) = In( sec (t)). 


LE mg0) cos (bsec(t) 1 
Au W(yi, 42) . 1 


Integrando temos que u>(t) = t. 
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e Uma solução particular é dada por 


Yp uy + Uzy> 


In( sec (?)) - cos (t) +t- sen (£) 


In( sec (In x)) -: cos (Inx) + Inx- sen (In x) 


Portanto a solução é 


y=c cos (Inx) + co sen (In x) + In(sec (In x)) - cos (Inx) +Inx- sen (In x) 


Exemplo 2.5.3. Considere a equação diferencial 


A-0)yl+xy -y=1-x 


(a) Verifique se y = e" e y> = x são soluções da equação homogênea correspondente. 


(b) Use o método de variação de parametros para determinar a solução geral da equação não 
homogênea. 
Solução: 


(a) y e y> são soluções da equação homogênea associada pois 


-d)yra,-n=U-nM(e)rate)-d =0 


A-dy +xy-y=(1-3)(0)+x(1)-x=0 
Portanto y e y> são duas soluções. 


(b) O Wronskiano de y e y> : 


Yi 42 e 
W= = =--xe=(1-ne 
, , X 
Ho e 1 
Para use o método de variação de parâmetros colocamos a equação na forma padrão 


dividendo por (1 — x) e temos 
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O método de variação de parâmetros fornece uma solução particular da forma y = 


u1y1 + U2y> onde 


jl= HG) sp fara] sRiE e 
DL WúyL,y) (I-de (A-ge = 


=X 





Integrando temos que m(x) = —-e “+ , dx (este integral não tem uma forma 





X mi 
fechada) 
u = st) = Ditos = 1 
2 uno) = Leg 
Integrando temos que u>(x) = — In(1 — x). 


Uma solução particular é dada por 


Yp = Uy + Uayo> 


= er . x — . 
— [-» + [Ea e Ia(l ==) 0% 


—x 
= a+e qd — xIn(1 — 2) 








Portanto a solução geral é 


=X. 


prod sex 1+e [ Eode-xinçi-a) 





Observação. O Exemplo 2.5.3 mostrou que em usando o método de variação de parâmetros podemos 
encontrar integrais muito difíceis de resolver. Além disso, determinar as soluções y e y> da equação 


homogênea com coeficientes variáveis também pode ser uma tarefa muito árdua. 


2.5.1.2 Método de Coeficientes a determinar 


O método de coeficientes a determinar é um método para determinar a solução particular 


yp da equação diferencial linear não-homogêneas de coeficientes constantes 


W + py +qu= (0) (2.5.9) 


quando f(x) tem a forma especial, envolvendo somente polinômios, exponenciais, senos e 


cossenos. 
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1º Caso: Quando a função g(x) = à) + ax + AO Edgar 


Admitiremos como solução particular a função 
Yplx) = Ag + Ajx + AR Psp dA 


Modificação: Se alguma parcela de y, coincide com termos de y então esse y, deverá ser 


multiplicado por x* onde s é o menor inteiro positivo que elimina essa coincidência. 
Exemplo 2.5.4. 

(a) Determine o solução geral da equação y” —-5y' +6=2x2—-1 

(b) Resolve o problema de valor inicial y' +2y = 122, y(0)=0, y(0)=1. 
Solução: 


(ay -5y +6=2x22-1 


Cálculo de yc : 
M=2>5y= 
Egn caraterística: A2-5A+6=0> e 
2 =35y=e* 
Portanto 
Ne= ce + coe* 
Cálculo de y, : 


Como g(x) = 2x2 — 1, assumimos a solução particular 
Yp =Axº+Bx+C. 


Como nenhum termo de y, coincide com termos de yc, não precisamos modificar yp. 


Para determinar os constantes 4, B,C, substituirmos y, na equação diferencial: 


YW =2Ax+B, y =2A 


Joseph N. A. Yartey & Simone S. Ribeiro 


114 


Substituindo na equação temos que 


2x2 —1 


2A—5(24x + B) + 6(Axº + Bx +) 


= (64)? + (6B — 104)x + (24 —5B + 60) 2x” = 1 


Comparando os coeficientes temos que 
1 
6A=2 3A= 3 


104 5 
Eine SE É 
FA RR 


MAs sDr6c==il o Db 2A +ob > 


6 ao 
Portanto 
55 
MES O dg O 
Solução geral é 
VV = Y+HYp 
2 
> RE do BDES 4 
= + + + +— 
c1º c2€ Po + 5 


(b)y” +2y = 12x 


Cálculo de yc : 

Ads gde 
Eqn caraterística: 122 +21=0> e 

A =25y=e* 
Portanto 

Ye=0+ Coe 
Cálculo de y, : 


Como g(x) = 12x2, assumimos a solução particular 


yp = AX +Bx+C. 
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Mas temos um termo constante em Yp e um termo constante em Yyç, logo modificamos 
yp pela multiplicamos por x assim eliminamos o termo constante. Portanto a expressão 
coreto de yp é 


Yp = AX +By+Cx 


Para determinar os constantes 4, B,C, substituirmos y, na equação diferencial: 
Y% = 34xº + 2Bx + C, Yo = 64x +2B 
Substituindo na equação temos que 


64x +2B+2(34%º +2Bx+C) = 12% 


= (64)xº + (64 + 4B)x + (2B +20) = 12x 
Comparando os coeficientes temos que 
64A=12 54A=2 


64+4B=0 B=-*E=-s 


2B+2C=0 5C=-B=3 


Portanto 


Yp = 2%é —-3x2 + 3x 


Solução geral é 


Yo = VerVo 


a+oe7+M IP +3x 
Agora usamos as condições inciais para determinar cy e c> : 


y(x) = cq + ce +23 +33 y(x) = —20e* +6% —6x+3 


Joseph N. A. Yartey & Simone S. Ribeiro 


116 


v0)=0 > q+o=0 
v(0)=1 > 20+3=1 


Resolvendo temos que cy = -1 e c, = 1. Portanto a solução do PVI é 


y="1+0 242 372 +3x 


2º Caso: Quando a função g(x) = Ea 


Admitiremos como solução particular a função 


Yp(x) = Ad 


Modificação: Se y, é uma das soluções de y. então esse y, deverá ser multiplicado por x onde 


s é o menor inteiro positivo que elimina essa coincidência. 
Exemplo 2.5.5. 

(a) Determine a solução geral da equação y” — 7y' + 12y =3e* 

(b) Determine a solução geral da equação y” — 7y + 10y = 8e” 

(c) Resolve o problema de valor inicial y' -4y +4y=8e€2, yO)=1, y(0)=5 
Solução: 


(a) 4-7 + Dy =3€* 


Cálculo de yc : 
M=35y= e 
Egn caraterística: 22-71 +12=0> e 
A =45y=e* 
Portanto 
Ve ce + ce 
Cálculo de y, : 


Como g(x) = 3º *, assumimos a solução particular 


Yp = Ae. 
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Como yp é uma das parcelas de yc, não precisamos modificar yp. 


Para determinar o constante A substituirmos y, na equação diferencial: 
yp=-Ae*, yp =Ae* 


Substituindo na equação temos que 


Ae*+7Ae*“ +1240* = 3* 


=> 2040” = 3* 


> A= ea 
20 
Portanto 
3 =X 
Yp = 50º 
Solução geral é 
yY = YetyY 


3x 4x 3 =X 
= qe +ce" + —e 
20 


(b)y” —7y + 10y = 8 


Cálculo de yc : 
A=2>n= 
Egn caraterística: 22 -7A+10=0> e 
Aa =5> p= 
Portanto 
Ve = cre + ce” 
Cálculo de y, : 


Como g(x) = 8e*”, assumimos a solução particular 
Yo = Ae, 


Mas temos que y, coincide com uma das parcelas de yc, logo modificamos y, pela mul- 


tiplicamos por x assim eliminamos este coincidência. Portanto a expressão coreto de yp 
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D+ 


= Ame 


Para determinar os constantes A substituirmos y, na equação diferencial: 
Y= Ae + 2Axe*, Y = 44€* + 4Axe* 


Substituindo na equação temos que 


44€* + 4Axe* — 7Ae* — 14Axe” + 104xe” = 8e% 
> 340 = 8e* 
8 
DA = —= 
3 
Portanto 
Yp = pes 
Solução geral é 
Vo = Y+tY 


8 
= qe" +0e* +20 — que? 


() y” — 4y' + 4y = 8e* 


Cálculo de yc : 
Apel pe” 
Egn caraterística: 22-41 +4=05 e 
A2 =25y =xe* 
Portanto 
Ve = cre? + coxe* 
Cálculo de y, : 


Como g(x) = 8e*”, assumimos a solução particular 
is AE, 


Mas temos que y, coincide com uma das parcelas de yc, logo modificamos y, pela multi- 
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plicamos por x e temos 


Yo = Axe 


Observamos que o modificado y, ainda coincide com uma das parcelas de yc, então 
multiplicamos de novo por x e assim eliminamos este coincidência. Portanto a expressão 
coreto de Yp é 


WE Ae 


Para determinar os constantes A substituirmos y, na equação diferencial: 
Y= 2Axe* + 24320, Y = 240” + 8Axe” + 8Ax2e” 


Substituindo na equação temos que 


24 + 8Axe + 4402 — 4(2Axe* + 24720) + 44x = 8% 
= 240 = 8 
>A =4 
Portanto 
w= Ads 


Solução geral é 


yY = Ye+Y 


= ce + coxe? + 4x2 
Agora usamos as condições inciais para determinar cj e c> : 


y(x) = (cy + cox + 422 )e2* = y'(x) = (c2 + 8x)” +2(c) + cox + 422 )e2* 


y(O) = 1 > cq=1 
Vad) =D E» e ES 


Resolvendo temos que c; = 1e c, = 3. Portanto a solução do PVI é 


v=(1+3x +42) 
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3º caso: Quando a função g(x) = sen (kx) ou cos (kx). 


Admitiremos como solução particular a função 


Yplx) = A sen (kx) + B cos (kx) 


Modificação: Se alguma parcela de y, coincide com termos de y então esse y deverá ser 


multiplicado por xº onde s é o menor inteiro positivo que elimina essa coincidência. 


Exemplo 2.5.6. Resolver as seguintes equações: 


(a)y” — 4y' +3y =3 sen (2x) 


(b)y"+y=4senx 


Solução: 


(a)y” -4y' +3 =3 sen (2x) 


Cálculo de yc : 
M=1>y=e 
Egn caraterística: 12-41 +3=0> e 
A =35y=e* 
Portanto 
Ve = e + Coe 
Cálculo de y, : 


Como g(x) = 3 sen (2x), assumimos a solução particular 
Yp = A sen 2x + Bcos 2x. 


Como nenhum termo de y, coincide com termos de yc, não precisamos modificar y,. 


Para determinar os constantes 4, B substituirmos y, na equação diferencial: 


y =24cos2x-2Bsen2x, y, =—4A sen 2x — 4B cos 2x 
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Substituindo na equação temos que 


(-44 sen 2x —4B cos 2x) —4(24 cos 2x —2B sen2x) = 3sen2x 


> (-A+8B)sen2x+(-B-84)cos2x = 3sen2x 


Comparando os coeficientes temos que 


-A+8B = 3 E 24 
5 A-= = eB=— 
B+84 = 0 bo bo 
Portanto 
3 24 
pec ea 2 25 “08 2x 
Solução geral é 
yY = YetYp 


= qe +ce*— = sen 2x — 8 cos 2x) 


(b)y'+y=4senx 


Cálculo de yc : 
A =i 
as 5 W = Cosx 
Egn caraterística: 1 +1=0> e > 
yo = senx 
Ao =—1 
Portanto 


Yo =C| cosx+co sen x 


Cálculo de y, : 


Como g(x) = 4 sen x, assumimos a solução particular 
Yp = Asenx+Bcos x. 


Mas temos o termo sen x está em y, e em yc, logo modificamos y, pela multiplicamos 


por x assim eliminamos esta coincidência Portanto a expressão coreto de yp é 


Yp = Ax sen x + Bx cos x 
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Para determinar os constantes 4, B, substituirmos y, na equação diferencial: 
Y = Asenx+Bcosx+Axcosx-—Bxsen x 


Y =2A cos x-—2B sen x— Axsen x — Bxcos x 


Substituindo na equação e simplificando temos que 
24cosx-—-2Bsenx = 4senx 


Comparando os coeficientes temos que 


Portanto 


Yp = —2X cos x 
Solução geral é 


yY = YetY 


= C| Cosx+c)senx-—2xcos x 


4º caso: Quando a função g(x) = g1(x) + --- gn(x) sendo que g;(x) é um polinômio (1º 


caso), ou exponencial (2º caso) ou seno ou cosseno (3º caso). 


Admitiremos como solução particular 


YplX) = Ypy +" Ypm 


onde y,, é a solução associado a g; respectivamente. 

Exemplo 2.5.7. Resolver as seguintes equações: 
(a)y” -9y = 2 + 5e* 
(b)y'-y=3senx-eZ+1 


Solução: 
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(a) y” —9y = + 5 


Cálculo de yc : 
MESES e 
Egn caraterística: 122 -9=0> e 
A =35yWy=e* 
Portanto 
Ve = 10% + coe * 
Cálculo de y, : 


Como g(x) = e2 + 5e*, a solução particular é da forma: 


Yp = Yp + Yp 


onde 
Yp, é uma solução particular de y' -9y =” e 


Yp, é uma solução particular de y” — 9y = 5e*. 


; 1 5 
Do 3º caso acima, temos que y,, = di e Y, = quer. 


Portanto 
1 ee > 3x 
=--€" + —xe 
Yp 5 6 


Solução geral é 


Y = Ye+Y 


a 1 5 
= qe*+ce — Ee + quer 


(b)y'-y=3senx-+1 


Cálculo de yc : 
Mel une 
Egn caraterística: 122 -1=0> e 
MA=-15y=e* 
Portanto 
Ve = ce +ce* 
Cálculo de y, : 
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Como g(x) = 3 sen x — e* + 1,a solução particular é da forma 


Yp = Ypy + Ypo + Yps 


onde 


Yp, é uma solução particular de y” - y=3 sen x, 


Yp, é uma solução particular dey' -y=- e 


Yp; é uma solução particular de y” -y=1 


Dos casos 1º, 2º e 3º acima temos que 
3 1» 
Yp = a senx, VYp = eae É” Yps = a! 


Portanto 


3 1 2x 
= "2 Eai —1 
Yp sen x e 


Solução geral é 


yY = VetY 


e 3 1 
= qe +ce *- sena- 5-1 


5º caso: Quando a função g(x) = g1(x) - g2(x) --: gn(x) sendo que g;(x) é um polinômio 


(1º caso), ou exponencial (2º caso) ou seno ou cosseno (3º caso). 


Construímos uma solução particular yp(x) da seguinte forma: 


(a) Escreve as soluções particulares y,; associada a cada função g; sem atribuir qualquer 


coeficiente, 
(b) Multiplique todos os y,, em (a) para obter todos os termos distintos do produto, 


(c) Insere os coeficientes a determinar na expressão na expressão em (b), um para cada termo, 


assim obtemos a forma correta de yp(x). 


(d) Caso um dos termos de yp(x) seja solução da equação homogênea associada, multiplique 


por xº onde s é o menor inteiro positivo que elimina essa coincidência. 


Exemplo 2.5.8. Resolver as seguintes equações: 
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(a) y” + 9y = (1º +27) 
(b)y"-6y +9y=e sen x 
()y' -2y — 3y = 2e* cos (32) 


Solução: 


(a) y” +9y = (12 +2)e* 


Cálculo de yc : 
Ay =3i 
= cos 3x 
Eqn caraterística: 122 +9=0> a as n 
y> = sen 3x 
Ao =-3i 
Portanto 
Vc = C1 COS 3X + 09 Sen 3x 
Cálculo de y, : 


Como g(x) = (x2 +2)e*” = g1(x) - g2(x) construímos a solução particular da seguinte forma: 


e Para gi(x) = Ds Ypi(X) = x +x+1 (sem os coeficientes) 


3x 3x 


e Para go(x) = e? — yp(x) = e* (sem os coeficientes) 
e Multiplicamos y,, € Yp, : 


(2 + x + 1 = 22% 4 xe + e 


Temos 3 termos distintos, assim teremos 3 constantes em Yp- 


Logo 
yp = (A + Bx + Cx) 


Para determinar os constantes 4, B, C substituirmos y, na equação diferencial: 


Yp = (A+2Cx +34 +3Bx +3Cx))e”, — yy =(2BC+3B+12Cx+3B+9A +9Bx +9Cx7)e* 
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Substituindo na equação temos os sistema 


C+3B+94=1 
19 1 1 
= A=— =-— = 
2C+3B =0 = 162” B 27" E 18 
18C=1 


Portanto 
fo fa Tas 
=|[— ——X+4— 
Jp (152 27" 3)e 


Solução geral é 


VV = Yc+HYp 
= Ccos3x+c sen 3x + ( a x + E oe 
= 8a e 162 27 18 
(b)y” —-6y +9y=e! sen x 
Cálculo de yc : 
W= Ex 
Egn caraterística: A2-61+9=05 ( AM =A,=3 5 
EA DADÃ 
y2 = xe 
Portanto 
Ve = cre + coxe 
Cálculo de y, : 


Como g(x) = e" sen x = g1(x) - g>(x) construímos a solução particular da seguinte forma: 


e Para gi(x) =" — yp(x) = é* (sem os coeficientes) 
e Para g2(x) = senx — Yyp(x) = senx+ cos x (sem os coeficientes) 


e Multiplicamos y,, e Yp, : 
e(senx+ cosx)=e senx+e cos x 
Temos 2 termos distintos, assim teremos 2 constantes em Yp- 


Logo 


yp = (A sen x + Bcos x) 
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Para determinar os constantes 4, B substituirmos y, na equação diferencial: 
y=e(Asenx+Bcosx+Acosx-Bsenx), y =€e'(24 cos x — 2B sen x) 


Substituindo na equação temos os sistema 


-4A+3B =0 3 4 
> a =-=—, B=— 
34+4B=1 25 25 
Portanto 
= [= sen x + E cos de 
Ip = 5 25 
Solução geral é 
VV = Yc+Y 
= (qto +(x sen x + = cos ae 
A 25 25 


Exemplo 2.5.9. Determine a função complementar y.(x) e uma forma adequada para a solução 
particular yp(x) utilizando o método dos coeficientes a determinar das seguintes equações 


diferenciais. Não avalie as constantes em yp(x) 
(a) y! -2y —3y = 2e* cos (3x) 
(b) y" + 25y = 4% sen (5x) — 2€* cos (5x) 
Solução: 


(a) y' —-2y — 3y = 2e* cos (32) 


AM =-1 a 
Eqn caraterística: A2-2A-3=05 5 Es Vi = Fi 
1=8 Ep 
Portanto 
Ve =c1€* + coe* 
Cálculo de y, : 


Como g(x) = 2e* cos (3x) = gi(x) - g>(x) - g3(x) construímos a solução particular da 
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seguinte forma: 


e Para gi(x) = É — Yp(X) = 2 +x2+x+1 (sem os coeficientes) 


e Para go() =" — Ypa(X) = e* (sem os coeficientes) 

e Para $3(x) = cos (3x) — Yyps(x) = cos (3x) + sen (3x) (sem os coeficientes) 

e Multiplicamos y,,, Yp, € Yp; : 
(+22 +x+ De cos (3x) + sen (35) = 9ºe* cos (3x7) +x2e* cos (3x) + xe” cos (52) + 
+e* cos (32) + 22º” sen (3x) + 22º” sen (3x) + xe” sen (5x) + e” sen (3x) 


Temos 8 termos distintos, assim teremos 8 constantes em Yp- 


Logo 


Yy = Age cos (3x) + Bx2e” cos (3x) + Cxe” cos (5x) + De” cos (3x) + Ex2e” sen (3x) + 


+ Fxe* sen (3x) + Gxe” sen (5x) + He” sen (3x) 


(Ax? + Bx? + Cx + De” cos (3x) + (Exê + Fx? + Gx + He” sen (3x) 


(b) y” + 25y = 4% sen (5x) — 2e* cos (5x) 


Cálculo de yc : 
MM =5i 
= cos (5x 
Egn caraterística: 2+25=0> e > e E 
yo = sen (5x) 
Ao = —5i 
Portanto 
Ye = C1 COS (5x) + co sen (5x) 
Cálculo de y, : 


Como g(x) = 4x? sen (5x) — 2e* cos (5x), a solução particular é da forma 


Yp = Yp + Yp 


onde 
Yp, é uma solução particular de y” + 25y = 4x sen (5x), 


Yp, é uma solução particular de y” + 25y = —2€* cos (5x) 
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Seja g1(x) = 4% sen (5x) um produto de um polinômio de grau 3 com a função seno. 
Construímos a solução particular y,, da seguinte forma: 

e Para 4º temos a solução particular 2 +x]+x+1 (sem os coeficientes) 

e Para sen (5x) temos a solução particular cos (5x) + sen (5x) (sem os coeficientes) 


e Multiplicando temos 


3 


(É + ea v( cos (5x) + sen (52) = 3º cos (5x) + x? cos (5x) +x cos (5x) + cos (5x) + 


+ sen (5x) + x? sen (5x) + x sen (5x) + sen (5x) 


Temos 8 termos distintos, assim teremos 8 constantes em y,, . 


Logo 


Yyp = AX cos (3x) + Bx? cos (5x) + Cx cos (5x) + D cos (5x) + Ex? sen (5x) + Fx? sen (5x) + 


+ Gxsen (5x) + H sen (5x) 


(Axº + Bxº + Cx + D) cos (5x) + (Ex + Ex? + Gx + H) sen (5x) 


Mas os termos D cos (5x) e H sen (5x) são soluções de y., então multiplicamos por x para 


eliminar esta coincidência. Portanto correta de y,, é 


Yp, = x(Azxê + Bxº + Cx + D) cos (5x) + x(Exº + Fx? + Gx + H) sen (5x) 
= (Axº + Bxê + Cx + Dx) cos (5x) + (Exº + Fxº + Gxº + Hx) sen (52) 
Seja g2(x) = —2e* cos (5x) um produto de uma função exponencial com uma função 


cosseno. Construímos a solução particular y,, da seguinte forma: 


e Para 2º temos a solução particular e* (sem os coeficientes) 
e Para cos (5x) temos a solução particular cos (5x) + sen (5x) (sem os coeficientes) 


e Multiplicando temos 


(> cos (5x) + sen (52) = e cos (5x) + e sen (5x) 


Temos 2 termos distintos, assim teremos 2 constantes em Ypo 
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Logo 


Ypo = Ie cos (5x) + Je” sen (5%) 


Portanto, 


= 
| 


Ypy + Yp 
= (AX! + By + Cx + Dx) cos (5%) + (Ex! + Fx? + Gx? + Hx) sen (52) + 


+ Te” cos (5x) + Je sen (5x) 


Resumo do Método dos Coeficientes a determinar: 
Para resolver y” + p()y' + q(x)y = g(x) 


(1) Encontre o conjunto fundamental de soluções [y1, >) da equação homogênea associada. 


Portanto a solução complementar é Wyc = ciy1j + Coy> 
(2) Divide, se necessário, g(x), em partes: g(x) = g1(x) + go(x) + :-- + gu(x) 


(3) Para cada g;(x), escolhe a sua forma da solução particular correspondente y,, de acordo 


com a tabela abaixo: 


( 
sen Bx ou cos px x(A sen Bx + B cos Bx) 


Pu(x)e”* sen Bx ou Py(x)e”* cos Bx (An ++ Age” sen px + (Bax” +---+ Bo)e”* cos Ex] 


onde s é o menor inteiro não negativo (s = 0,1,2,...) que garanta que nenhuma parcela 





da solução particular seja solução da equação homogênea associada. 


(SD Yp =Yp + Yp Ho + Yp 
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2.5.1.3 Exercícios 


1. Resolva a equação diferencial dada pelo método dos coeficientes a determinar. 


(a) y' +3y +2y=6. Resp:y=ce 2 +ce*+3 

3 
(b) y” -10y+25y=30x+3. Resp: y=(ci+ ce” + = (2x + 1) 
(c) av +y +y=2-2x. Resp:y=(g+oxne?+x-4x+5 


-2 8x? 
(d) y' +3y = -A8x2e*. Resp: y=ci+ce + = + Ee 


(o) yw'-y=-3 Resp:ry=g+tod+3x+3 


O y-y+qu=3+ e? 


1 
- Respy=(c + cone +3 + pve 
3 
(g) y' +4y=3sen2x. Resp: y=c4 cos (2x) +» sen (2x) — 7* cos (2x) 


1 1 
(h) y'+y=2xsenx. Resp: y=cicosx+c> senx+ x senx— x cos x 


2. Resolva a equação diferencial dada sujeita às condições iniciais indicadas: 


(a) y'-6y —-7y=-9*, y0)=-2, y(0)=13. Respiy=e”-2" -e* 
(b) 5y" +y =-6x, y(0)=0, y(0)=1 


(o) y'+4y +5y=35€"*, yO0)=-3, y(0)=1 
2 
(d) a +alx=Fosenwt, x(0)=x(0)=0 


(e) y' +9y= cosx— senx, y(n/2)=y(m/2) =0 


3. Resolva a equação diferencial dada pelo método da variação dos parâmetros. 


(a) y" +y= secx 
(b) y"+y= senx 
(0) y' +y= cos?x 
(d) y” —-y=coshx 
(e) y' +3y +2y = Ro 
(DD y'+2y +y=e"lInx 


4. Resolva a equação diferencial pelo método da variação dos parâmetros, sujeita à condição 


inicial y(0) = 1, y'(0) = 0. 
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(a) 4y” E y — xe*/2 


(b) y' +2y —-8y=20* -e* 
(0) 2y" +y —-y=x+1 


(d) y' — 4y' + 4y = (12% — 6x) 


5. Determine a função complementar yc(x) e uma forma adequada para a solução parti- 
cular yp(x) utilizando o método dos coeficientes a determinar das seguintes equações 


diferenciais. Não avalie as constantes em y,(X) 


33% — 26X sen 9x 


Ve = cre + coe* 


yp(x) = (Axº + Bx! + Cx + Dx? + Ex) + Fe cos 9x + Ge * sen 9x 


(a) y” —-9y=6x 


Resp: 


(b) y” -4y' +4y=52 -2-— 2202 + 4€2 cos x 


Ve = cj + coxo” 


Resp: 
yplx) = AÉ +BX +Cx+D+(Ext+ FÃ + Ge + He” cosx+ Ie sen x 
(O) y' +4y +20y=xe * sen4x -3 cos 4x — xe 2 
Vc = ec cos 4x + co sen 4x) 


Resp: 4 yp(x) = (As? + Bxº + Caje * cos 4x + (Dx + Ex? + Ex)e * sen 4x + G cos 4x+ 
+H sen 4x + (Ix + De” 


1/2 1/2 


6. Sabendo-se que as funções “sent e t”"'“ cost são soluções linearmente independentes 


1 
da equação Px+tr+ (? — 5) x=0, t>0,encontre a solução geral de 
1 
Pr+tr+ (2 E 5) = 382 sen t. 


7. Determine duas soluções linearmente independentes de 2x — 2x = O da forma x = f”. 


Usando essas duas soluções, determine a solução geral de Pr-2x=P. 


8. Uma solução da equação x + p(t)x + q(t)x = 0 é (1 + t2, e o Wronskiano de duas soluções 


qualquer, desta equação, é constante. Determine a solução geral de 


x+p(x+ag(tx=1+t 


MATB97 - Equações Diferenciais 


133 


Aula 2.6 


Equações diferenciais lineares de ordem 


superiores 


Agora vamos começar a construção de soluções para equações diferenciais lineares de ordem 


n, isto é equações diferenciais da forma: 


d"y dry dy n 
MO ar + Ma) esse a + ag(x)y = g(x) (2.6.1) 


A equação é 
e de coeficientes constantes: quando po,p1,:-: , pn são todos constantes 
e homogênea: quando g(x) = 0 
e não-homogênea: quando g(x) + O 


e homogênea: quando g(x) = O 


forma padrão: divide (2.6.1) por a,(x) 


d"y dy dy o. 
dr + Pi a + + Pao) o + pnlody = 8000) (2.6.2) 





Definição 2.6.1. Uma solução de ( 2.6.2) é uma função y = h(x) definido no intervaloI c R que 
tem derivadas y =h'(1), y' =h"(9),-,y” = (x) e satisfaz a equação (2.6.2) para todo x 
no intervalo 1. 


tr 


Exemplo 2.6.1. Mostre que xe * é uma solução da equação y“ +2y” +y” = 0 paratodox e R. 


Joseph N. A. Yartey & Simone S. Ribeiro 


134 


Solução: 
use 
então 
= cg t ==” 
preto je =D er 
Up te Dep = See 
WO=-e*-(3-me*=(-4+me* 
Portanto 


VS ty tp Ad One" it met= (Ars +r6=21-2 +30" =D 


Portanto y = xe * é uma solução. 


Problema de Valor Inicial (PVI) 
Um problema de valor inicial da equação diferencial linear de ordem n consiste na equação 


d"y dra, dy n 
qa pi) am prior Pro + pnlc)y = g(x) 


en condições iniciais 





Y(xo) = yo 

dy ; 

da (0) =Y 

dy (n-1) 
dyxn-1 (xo) MO 5 
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Teorema 2.6.1 (Existência e Unicidade). Se as funções pi(x), p>(x),::-,pnlx) e g(x) são 


continuas e deriváveis no intervalo aberto (a, b), então existe somente uma função y(x) satisfazendo 


(2.6.2) no intervalo (a, b) e as condições iniciais 


Y(xo) = Yo 





Soluções da equação linear homogênea de ordem n 


Definição 2.6.2. Um conjunto de funções lp, P>,:-- Pn) é dito linearmente independentes no 


intervalo 1 = (a, b) se não existem constantes c1,c>,:-* ,Cn tal que 


cipi(x) + copa(x) +: + cnPnlx) = 0, Yxel 


excetocy =c)=:::=c,=0. 


Teorema 2.6.2. Um conjunto de funções deriváveis ly, P>,::* Pn) é linearmente independentes 


no intervalo T se e somente seo Wronskiano do conjunto é diferente de zero, isto é 


0 * Wíg, Pa, od » Pn)(x) = det 
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Teorema 2.6.3. Suponha que as funções pi(x), po(x),-:- , pn(x) e g(x) são continuas e deriváveis 


no intervalo aberto (a, b) e as funções yi(x), y>(x),:-- , YnlX) são soluções da equação homogênea 


d"y 
dx” 


dr ty 


dyxr-1 


+ pi(x) 


d 
Ea de pao + pu(x)y =0. (2.6.3) 


Então se W(y1, y>,:** , YUn)(x) £ O em pelo menos um ponto a < x < b, então a solução geral de 
(2.6.3) é dado por 


Ve = Caya(x) + coyolx) +: + CnYyn(x) 


c1,C2,º:* ,Cn constantes arbitrárias. 





Soluções da equação linear não-homogênea de ordem n 


Considere a equação diferencial linear de ordem n não-homogênea 


d"y 
dx” 


diria; 
dxn-1 


d 
eae prod + pa(o)y = 8(2) (2.6.2) 








+ pi(x) 


e suponha que y(x), y>(x),--* , Yn(x) é um conjunto de n soluções linearmente independentes 
da equação homogênea correspondente. Se yp(x) é qualquer solução particular de (2.6.2) então 


a solução geral de ( 2.6.2) é da forma 


VC) = yp(0) + Velo) = YplX) + cryi(a) + coyalx) +: + Cnynlx) 


2.6.1 Equações Diferenciais Lineares Homogêneas de ordem n com 
coeficientes constantes 


Consideramos a equação diferencial linear homogênea de ordem n com coeficientes cons- 
tantes 
d"y dry 


dy 
ni O a o dO (2.6.4) 


Quando n = 2 sabemos que uma solução desta equação pode ser resolvido, assumindo inicial- 


mente uma solução da forma 


u(x) = e!* (2.6.5) 
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e substituindo e resolvendo por À. Vamos fazer a mesma coisa para o caso geral n. Substituindo 


(2.6.5) em (2.6.4) obtemos o polinômio característico: 
A" + Rear Sfpstoms AA? +01A+09=0 (2.6.6) 


As raízes do polinômio característico determinam a solução de ( 2.6.4). Em relação o polinômio 
característico, tem-se 3 casos a considerar: 
Caso 1: O polinômio característico admite somente raízes reais e dis- 
tintas 

Sejam A1, À>,:-- , An as raízes do polinômio caraterístico, então temos n soluções linearmente 


independentes 


À An-1X AnX 
r 


nn) =, plo) =", e palj)=e Yn(x) = 


Assim a solução geral é 


v= ese +. + cy en + cel 


Exemplo 2.6.2. 
(a) Determine a solução geral da equação y 9 — 13y” + 36y = 0. 


(b) Resolve o problema de valor inicial y” -5y' -22y +56y=0, y(O)=1, y(0) =-2, (0) = 
—4, 


Solução: 


(a) Polinômio caraterístico: A4-1322+36=0 


2 =951=3,A=-3 
2 13+ V169-144 13+5 | 


A 2 2 


ou 


2 =4514=2,M=-2 


Portanto solução geral é 


v= cet + ce + ce? + cye 
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(b) Polinômio caraterístico: p(A) = A -5A2-221+56=0 
Primeiramente observamos que p(2) = 8 — 20 — 44 + 56 = 0 portanto À — 2 é um fator de 


p(A). Aplicando o Briot-Ruffini 


1 -5 —-2 56 2 





1-3 -28] 0 
Isto nos dá a fatoração 


p(A) = A2-5A2-221 +56 
= (A-2(A?-3A- 28) 


Mas A?-31-28=(A-T(A+4) 


Portanto as raízes de p(A) são À =2, A) =7, À3 = —4. Portanto a solução geral 
y= ce + co + cao * 


Agora usamos as condições inciais para determinar c1, ce c3: 


Vos ce + coe”* + cao 
=> y(x) = 201€* + 7c,e* — 4c3e “E 
> y'(x) = 42 + 49c2e” + 16c30 ** 
y(O) = 1 > cq+o+c=1 


v(0)=-2 > 24 +70, — 403 = —2 
v(0)=-4 > 40 +49, + 1603 = —4 


Resolvendo temos que | = —, o =-— e c3 = —. Portanto a solução do PVI é 


13 16 14 
is 55 33: 


13 2x 16 7x 14 —4x 
E / vim — — +— 
1 550 295" 
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Caso 2: O polinômio característico admite raízes reais com alguns 


raízes com multiplicidade. 


Suponhamos que o polinômio característico possui uma raiz real A, com multiplicidade 
p. Com essa raiz construiremos p soluções linearmente independentes da equação (2.6.4) da 


seguinte maneira: 


=: AX ES as E ARS ÇÃO — p-2,AkX =p ,ÃRX 
y=e*, y=x*, yy=xe 8, ypy = yy = cet 
Essa construção é feita com todas as raízes reais do polinômio característica, para obtemos 


n soluções linearmente independentes da equação ( 2.6.4). 


— ApaX — Ap42X = 6AnAX — AnX 
Yp+i = err, Yp+2 = er, “e, Ynm = err, Yn=e" 
Assim a solução geral é 
v = gts cxd + carte po. + Cpo e + Ene tea + Eterno + Epa ar ++ cnetrt 


= (g+ox+cr++ cp ete + Copgelreit ++ cnetrã 


Exemplo 2.6.3. 

(a) Determine a solução geral da equação y 9 + y” -3y' -5y -2y=0. 

(b) Resolve o problema de valor inicial y) -y' -y +y=0, y0)=2, y(0)=1, (0) = 0. 
Solução: 


(a) Polinômio caraterístico: p(A) = M+A2-342-5A-2=0 
Podemos verificar que 
p(-)=1-1-3+5-2=0 
p2)=16+8-12-10-2=0 


Portanto —1 e 2 são raízes de p(1), podemos aplicar Briot-Ruffini 
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Isto nos dá a fatoração 


EA == Bj=? 


ni 
Rs 
I 


(A+I(A-2(A2 +21+1) 
(A+ IXA=2AA + 1) 


(ADA —2) 


Portanto as raízes de p(A) são À; = —1, (multiplicidade 3) e À, = 2. 
Portanto a solução geral 


v=e%+(c+xc3 +12)” 


(b) Polinômio caraterístico: p(A) = A -A2-A+1=0 
Primeiramente observamos que p(1) =1-1-1+1 = 0 portanto À — 1 é um fator de p(1). 


Aplicando o Briot-Ruffini 





Isto nos dá a fatoração 


flo) = AÉ-22-A+41 
= (A-1)(42-1) 
= (A-I(A-I(A+1) 
= (A-I(A+1) 


Portanto as raízes de p(A) são À4 = 1 (multiplicidade 2)e A) = 1. 
Portanto a solução geral 


y=(ci+xc)e + cge* 


Agora usamos as condições inciais para determinar c1, ce cs: 


u(x) = (g+xc)e +cge* 
> y(s) = (c+co+xco)e — ce” 
> y'(x) = Qo+ca+xco)e +cae* 
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y(O) = 2 > q+c=2 
VO=1 = qro gel 


v(0)=0 > 29+c+c3=0 


Resolvendo temos que c; =2, c) = -1 e c3 = 0. Portanto a solução do PVI é 


y=(2- ne 


Caso 3: O polinômio característico admite raízes complexas distintas. 


Observe que as raízes complexas aparecem em pares conjugados. 
Para cada par distintas raízes complexas conjugados (À, À), À = a + ip (com multiplicidade k) 


podemos construir 2k soluções linearmente independentes da seguinte forma: 


wm=e"cos(pa), y=xe“cos(pi), y=2e“cos(pro, , yp= Te cos (px) 


visi = € sen (Bx), Y+2 = xe” sen (Bx), Yws = x2e“* sen (Bx), cc, Yk= 216% sen (Bx) 


Exemplo 2.6.4. Resolver as seguintes equações 
()y9 -y=0 
(b) y + 12% + 1049) + 408y" + 1156y =0 
Solução: 


(a) Polinômio caraterístico: p(A) = A*-I=(A-DA+-DA-A+LD=0 
Portanto as raízes de p(A) são A = 1, À =—1, Ag =1, As =. 


E podemos construir as seguintes soluções linearmente independentes: 


A =1 > y=e 
Ao =-1 > Yy= e” 
== ; va(x) =P cos x = cos x 


va(x) = e" senx = senx 


Portanto a solução geral 


y=ce+ce*+cycosx+cysenx 
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(b) Polinômio caraterístico:  p(A) = Àº + 121º + 10412 + 4081? + 11561 = 0 
Mas 


Aº + 121º + 10422 + 408122 + 11561 = A(A!+ 122º + 1041? + 4084 + 1156) 


A(A? +61 +34? 


= MA+3+5)(A+3- 5? 
Portanto as raízes de p(A) são 
A =0, A2=-3-5i, (multiplicidade 2), A, = 3 + 5i (multiplicidade 2) 


E podemos construir as seguintes soluções linearmente independentes: 


M=0 > p=e&=1 
va(x) = e * cos 5x 
= va(x) = xe cos 5x 
Ao =-3-—51, À = —-3+ 5 (multiplicidade 2) — 
va(x) = e sen 5x 
ys(x) = xe ** sen 5x 


Portanto a solução geral 


y=ci+ce* cos (5x) + ce ** sen (5x) + cgxe * cos (5x) + coxe ** sen (5x) 


2.6.2 Equações Diferenciais Lineares Não-Homogêneas de ordem n 
com coeficientes constantes - Método dos Coeficientes a deter- 
minar 


A solução geral da equação diferencial linear não-homogêneas de ordem n com coeficientes 


constantes é 
dºy dy dy , 
o + MATE ER ax + aoy = g(x). (2:67) 
é dado por 
YO) = Yp + Vc 
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onde yc é a solução geral da equação homogênea associada e y, é uma solução particular de 


(2.6.7). 


A resolução de y, será pelo Método dos coeficientes a determinar que foi estudada na Aula 


12. 
Exemplo 2.6.5. Resolve a equação diferencial ou o problema do valor inicial 
(a)y” -4y =1-3x 
(b) y” —3y" +3y' =y = 10xe' 
(c) y” — 4y' = 12 sen 2x 
Solução: 
(a)y” -4y' =1-3x 
Cálculo de yc : 


A =0>5ym=1 


Ao =25 2 = e 
Egn caraterística: 2 -41=0> E 
e 


Ag=-25y= ga 


Portanto 


Ve=0+ co + cao ** 


Cálculo de y, : 


Como g(x) = 1 — 3x, assumimos a solução particular 
Yp = Ax+B. 


Como o termo de y, é uma solução em yc, precisamos modificar y,. pela multiplicação 


por x. Portanto a forma correta de y, é 
Yp = Ax? + Bx 
Para determinar os constantes 4, B substituirmos y, na equação diferencial: 


YW =24x+B, y=2A y =0 
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Substituindo na equação temos que 


H 
pa 
| 
Q 
Q 


(0) — 4(2Ax + B) 


> 84x —4B 


H 
pa 
| 
Q 
Q 


Comparando os coeficientes temos que 
-8A=3 3A= e 
5 
Spies 
o 4 


Portanto 


Solução geral é 


pe. E De Um 


a 3x? 
q+0e+ce + a 


REIS 


(b) ins —3y” +3y' E y — 10xe* 


Cálculo de yc : 
Egn caraterística: 42 -3122+31-1=0>(A-12 =0> A =1 (multiplicidade 3) 
Portanto 
Ve = ce + coxe' + care 
Cálculo de y, : 


Como g(x) = (10x)e* = g1(x) - g2(x) construímos a solução particular da seguinte forma: 


e Para gi(x) = 10x — yp(x) =x +1 (sem os coeficientes) 
e Para gx(x) =" — Yyp(x) = é” (sem os coeficientes) 


e Multiplicamos y,, e Yp, : 


(x+ e = xe" +e 


Temos 2 termos distintos, assim teremos 2 constantes em Yp- 
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Logo 


Yp = Axe" + Be” 


Mas os termos Axe" e Be" são soluções de y., então multiplicamos por xº para eliminar 


esta coincidência. Portanto correta de yp é 
Yp = 2 (Ax + BJ 


Para determinar os constantes 4,B substituirmos y, na equação diferencial temos os 


sistema 
244 = 10 5 
> (a =—, B=0 
B=0 te 
Portanto 
5 
Yp = o x* e 


Solução geral é 


= Ver Up 


2.6.3 Exercícios 


1. Ache a solução geral de cada uma das equações diferenciais abaixo: 





(a) y” +3y' -4y —12y=0. Resp:y=cet+ce+ce 
(b) 2y9 +11y” +18y" +4y -8y=0. Resp: y=ciel2 + (co +cax+ cg?) 
(c) vu +16y=0. Resp: y= ce + Coe” + c3 cos 2x + cy sen 2x. 
(d) y” —-5y"+7y' =0. Resp:y=c+[|c cos [3] +c3 sen [E] dz. 
(e) 49 — 1549 + 84y9 — 220" +275y — 125y=0. 
Resp: y= e" + ce + caxe* + cge” cos x + c5e”* sen x. 
(£) vt — 5y" +4y = E ra, Resp: y= ce! +ce”* + Es pg = no + res 


(MD) y”!-3y' +3y -y=e, Resp:yv=(g+0ox+cad)e + 
8) Y y Voy py 


(h) y” -2y =322 -2x+1. Respiy=c+ox+ce*-— — — — a 
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1 
(i) vt o 4y" = 8sen 4x. Resp: yV=0 + C9X + cao + cge?* + Z0 sen (42). 
() vt + 2y” + y= 4 sen X. Resp: yY = Cc Cos x + c) sen X + CgXx COS X + C4X Sen X — 
1 
2 sen (x). 
(6) 4) -5y" +6y =2x+3e2. Resp:. 


1) y)+4y =3* -x-2sen2x. Resp:. 


2. Ache a solução dos seguintes problemas de valor inicial: 


tr 


(a) y” +y = senx, y(O) = y'(0) = y“(0) = 0. 


Resp: y=-xcosx+3x-—2 senx 
(b) y -6y9 -y! =72x+24, Y(0) = 0,y(0) = 0, y(0) = 6, y(0) = —57. 
Resp: . 


3. Resolve as equações de Euler-Cauchy de ordem 3 


(a) xy” +32y" -2xy +2y=0 Resp: cix+cox +caxlnx 
2 
(b) xy” +2xy —2y =x Inx+3x Resp: cx + cx cos (In x) + c3x sen (In x) + 5 In x — 


x +3xlnx 


1 
() 2y” +xy -y=xmx Resp: cx+coxlnx+ cax(In x)? + q (In ae 


4. Determine a função complementar yc(x) e uma forma adequada para a solução parti- 
cular yp(x) utilizando o método dos coeficientes a determinar das seguintes equações 


diferenciais. Não avalie as constantes em y,(x) 


(a) 9 +2y" +y =x+xcosx 
Vc=C1+C2 COS X+C3 SeNnx+c4xCOos X+C5xsenx 
Resp: 
yp(x) = x(Ax + B) +32 [(Cx + D) cos x + (Ex + F) sen 2] 
(b) yP + 8y” + 16y = sen 2x +x2 cos 2x 
E Vc = C1 COS 2X + C2X COS 2X + C3 Sen 2x + cgx sen 2x 
esp: 
Yp(x) = (4º + By + + Di) sen 2x +(Exº + Ex" 4 Ga + Hx) cos 2x 
(0) VP +y' =x cosx 
Ve=C|+C2X+C3 COS X+C4 sen x 
Resp: 
Yp(x) = x(Ax? + Bx + C) sen x + (Dx? + Ex + F) cos x 


(d) y” +8y' + 16y=xsenx+x2 cos 2x 


(e)  déá —3y" +3y —-y= xe” —- 3 
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(8 VP +2y"+y=xsenx 


(g) A = 438) +6y' -4y +y= Pere 
Rss: Ve = e + coxe + cale + cute 
yplx) = x(A + Bx + Cx? + Di)e” + (Ex? + Fx + G)e* 
(h) yP +5y' +4y= senx+ cos2x+ sen3r 


(1) y”+2y"+2y =3"* cosx 


pass Ve=cC+ce"cosx+ce*senx 
e Yplx) = x(Ae* cos x + Be* sen x) 
Du! +2y"'+2y =xe*cosx-xe” senx 
(9) yP +4y%) + 8y" +8y +4y=7e* cos x 
Rd Vce=ce"cosx+cxe*cosx+c3e “senx+caxe* senx 
a yp(x) = x2(Ae”* cos x + Be” sen x) 


OD y9-2y' +y=xe + 20% + e 
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Aula 2.7 
Transformada de Laplace 


2.7.1 Definição e Exemplos 


Definição 2.7.1. Seja f : [0,+00] > R um função. A transformada de Laplace de f(t) denotado 
por Llf(t)) ou F(s) é definido pela integral 


09 A 
LI) = F(s) = ij ES FO dte= fim [o ef dt 


Ao 


desde que a integral convirja, onde s é um parâmetro (real ou complexo). 


Teorema 2.7.1. [Condições Suficientes para a existência de LIf(H] Se f : [0,00] > R é uma 


função continua por partes e existem constantesa >0, M>0€eT>0Otal que 


ICI] < Me” para todo t>T, 


então Llf(t)) existe para s > a. 





Exemplo 2.7.1. Determine a transformada de Laplace da função f(t) = t”, t>0, n= 


0,1230004 
Solução: 


e Paran =0, 


A -st |Á 
£t1) = lim et1dt= lim É | =lim A fes — 1] = o (s > 0). 


ASoo 0 AS5o —S 0 ASo S 
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e Paran>1 





A 
Lt”) = lim eo edi 
Ao Jg 
A 
pf. est A nti-le-st 
= lim se ii ——— dt (Integrando por partes) 
ASo0o —8 0 0 =o 





A 
Ns = é 
=> 0+- lim E pd 
S AsSo0o 0 


n 
= qt), (s > 0). 
Portanto temos a relação de recorrência, 
n n n—1 
Lf)= Lt) Yn, 


que significa 











= 1 E 2 
Lp! a Lo er, fps ao Er. o 
Potanto por indução, temos que que 
Ley = Eeury=0.0 À gn?) EA no? ps) = 
8 S s s s s 
= as n=1 n-2. ops E 
s 5 S s s s sn 


Portanto temos que 


cty= E (s>0. 


gn+1” 


00 A 
Observação. A partir de agora vamos usar / em vez de lim A , Sem prejuízo. 
0 Ao Jg 


Exemplo 2.7.2. Determine a transformada de Laplace de f(t) = €*, a ER. 
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Solução: 


OO 


e 4 et dt 


Le! = 


e(smadt dt 


Ss 


eTlsmat 


—Ls — a) 





t=0 


>, s>a. 
s-a 


Exemplo 2.7.3. Determine a transformada de Laplace de sen (at) e cos (at). 


Lembrando que 


AX E Er XxX — ES 
[Je sen (bx) dx = ar sen (bx) — b cos (bx)), [| e“ cos (bx) dx = ap cos (bx) + b sen () 











Solução: 
Li sen anj= [ es. sen (at) dt Li cos anj= [ es. cos (at) dt 
0 0 
=8E io et ço 
E (—s sen (at) — a cos (at)) = = ZrAt cos (at) + a sen (at)) , 
a s 
— RR s>0. = SR! s>0. 


Exemplo 2.7.4. Determine a transformada de Laplace de sinh(at) e cosh(at). 





Solução: 
Lisinh(at)) = N es. sinh(at) dt 
0 
E pr est o (et E e) E” 
0 2 
= | et. et a > [— Ed 
2 Jo 2 Jo 
1 1 1 
-s — = | (por Exemplo 2.7.2) 
a 
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8 


Llcosh(at)) es! . cosh(at) dt 


Il 


SS 


et + et 
est. e) 3 dt 





1 [9 1 [O 
= |, eseetars [ Es e di 
2 Jd 2d 
1[ 1 1 
>= |— + — E lo 2.7.2 
ls el ops ) 
o s 


Tabelando os resultados acima temos 


Fe = 21/60) 
E n! 
1 
at 
sen O 
Ss 
a 

Ss 

Sa 
Ss 

n2” 


s2 ey 
s2 — Cu 


Tabela 2.7.1: Transformada de alguns funções 





Lema 2.7.2 (Linearidade da Transformada de Laplace). Se a e 6 são constantes, então 


Lta ft) + Bet) = a fB) + BLIg(b) 


para todo s tal que as transformadas tanto de f quanto de g existam. 
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Demonstração. 


| ea f(t) + Ba(t) ját=a [ e Ft) arg e g(t) dt 


sempre que ambos as integrais convergirem para s > c. logo, segue que 


af) + BLig()) 
aF(s) + BG(s) 


La f(t) + Bgtt) 


Exemplo 2.7.5. Calcule 


(a) L(1 + 5t) (b) £L/4e*! — 10 sen 2) 





Solução: 
(a) LIL +58 = LH) +5LI) = -+ = pela Tabela 2.7.1 
O LUSA Ns dn sida sense =.  peatbaa do, 
s+3 8244 


2.7.2 Transformada Inversa de Laplace 


Teorema 2.7.3. Se f(t) e g(t) satisfazem as hipóteses do Teorema 2.7.1e 


Lig(b) paratodo s >a, 


g(t) em todos os pontos t onde ambas as funções são continuas. 





é única. Isto permite definir a Transformada Inversa de Laplace £”! como 
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Um método conveniente para se obter as transformadas inversas de Laplace, consiste em 


usar uma tabela de transformadas de Laplace. 


Exemplo 2.7.6. Da Tabela 2.7.1 temos que 


Lema 2.7.4 (Linearidade da Transformada Inversa de Laplace). Se a e B são constantes, então 


LHaF(s) + BG(s)) = «LL HF(s)) + BL HG(S) 


onde F e G são as transformadas das funções fe g. 





Exemplo 2.7.7. Calcule 


E Ê (o) 17 [ES] 


st s244 s2+4 








Solução: Usando a linearidade da transformada inversa de Laplace e Tabela 2.7.1 temos que 


(a) 17 (5 E p=) E |=P - sen (28) 


st 2 +4 st s2+4 
—2s + 6 —2s 6 —2s 6 
o ER) = ctfgiresta) = (gere (sb] = 2 sen 
(JL s2+4 L sS2+4 s244 E s2+4 L s2+4 copa) 
3 sen (21). 




















Transformada Inversa de Laplace e Frações Parciais 


Frações parciais desempenham um papel importante na determinação das transformadas 
inversas de Laplace. 

Se uma transformada F(s) não puder ser encontrada na tabela, então deve-se expandir em 
frações parciais e escrever F(s) em termos de funções simples de s nas quais as transformadas 


são conhecidas. 
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Exemplo 2.7.8. Determine a transformada inversa de Laplace das seguintes funções: 
1 3s + 1 
a —1 
——— b —— = — p. 
Gs (emos) E (e) 


Solução: 


(a) Usando frações parciais temos que 


É aÃ oB 
(s+1(s+2) s+1 s+2 





1 1 
Resolvendo temos que 4 = e B= ER Portanto 


Lo IB 18 


(s+1(s+2) s+1 s+2 








Logo, 


- 1 sap BY pa[18B)Y ha 
e ima! il E (53) 31º Pe, 


(b) Usando frações parciais temos que 


3s+1 A Bs+C 


—D—D— =" + ; 
(s-1s2+1) s-1 841 





Resolvendo temos que 4 =2, B= -2eC=1. Portanto 


3s+1 o 2 1 2s 


DD 2 +>— 
(s-1(s2+1) s-1 8241 58241 


Logo, 
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2.7.3 Propriedade da Transformada de Laplace e Sua Inversa - Pri- 


mero Teorema da Translação 


Lema 2.7.5 ( Primeiro Teorema da Translação). 


(a) Se Lif(t)) = F(s) ea for um número real qualquer, então, 
LH) =F(s—-a) ou Le FB) = LfBlsssa- 


(b) Se L!(F(s)] = f(t) então 


LHF(s = a) Er (£) ou EO) sasa = er ) 





Demonstração. 
(a) 


Len) = a et et ft) dt 


E Ki Cet ft) dt 
0 


| et f(t) dt, (pela mudança u=s-—a) 
0 


= Eu = FS = a): 


(b) Seja F(s) = [ e f(t) dt então 


0 


F(s-a) = E ú eo f(t) dt 


= e ee ft) dt 
0 
= Me ft). 


Portanto, L!F(s — a)) = e“ f(t). 
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E 
Exemplo 2.7.9. Calcule 
(a) Lte ?! sen 3) (b) Lt E) 
—1 1 —1 1 —1 2s + 5 —1 s + 1 
a tos fo, (25573) ia (E 0) ras E 
Solução: 
(a) Lle “* sen 3t) = Li sen 3t)sss-(-2) = 215 o Go 
531 PA me 3 
(b) Lle t J = Lt Jsss-(5) = 4 que = (5 5 
4 1 
o [e=3] 
A 


Seja F(s) = = então f(t) = LF(s)) = a 





Agora G — = F(s— 5). Portanto, 
—1 1 = £ (F(s o 3) Es e f(t) = 1 açã 
(s— 3) 410 
(dj [=] eo bp completando o quadrado 
2-45+5) (6-22+1 P q 


Seja F(s) = + então f(t) = L F(s)) = sent. 


= F(s — 2). Portanto, 


1 
A Ego 
EREs (s-22+1 


mm —+ = L!F(s-D= ft) = e sent. 





=| 2ssPo 
ler) 


Usando frações parciais, temos que 


2+5 2 1 


—— + 


(s+22 s+2 (s+22 
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Portanto 














do [Eds Do MIEgUiR SD E 1 E 1 
L (mL I=5)+£ om ns | 


Seja F(s) = = então f(t) = L!(F(s)) =t. 





Agora +) = F(s + 2). Portanto, 
-1 1 “pa = It 
£L o =L'(F(s+2)=e“f(t)=te*. 
Portanto 











= | ZE | am) Ro “1 1 = eb qi 1 E. 
e (EL (rat om) o om) RR 


(f) £ >= nl Es —— TZ TI6 completando o quadrado. 


Temos que 
L (ra) + (ES) 4 [sr] ci lesma! 


+3 
e Vamos calcular £! (Ss! 





Seja F(s) = 5 - 5» então f(t) = L!F(s)) = cos (4). 
+3 = 
Agora 6+32+42 = Fís E 3). Portanto, 


E EE | =1 (+39) =) = e cos (49) 


1 
e Vamos calcular £! rsss! ; 


sen (4t) 
A 





Seja G(s) = então g(t) = L! (G(s)) = 


1 
s2 + 42º 
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Agora E = G(s +93). Portanto, 
eg fotus Crer 3 sen (4) 
E sr] Sd CRASE se 
Portanto, 
=1 ca a - pd BED E = 1 ne. Edi sen (4t) 


2.7.4 Transformada de Laplace da Função degrau unitário 


Funções degrau unitário 





Transformada de Laplace da Função degrau unitário 


Exemplo 2.7.10. Determine a transformada de Laplace da função degrau unitário U(t — a). 
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Solução: 


Mu(t-a)) = a est.Ult-a)dt= es dt 
0 a 








Expressando funções definidas por mais de uma sentença em termos 


de funções degrau unitário 


Podemos expressar uma função f(t) definida por mais de uma sentença como uma combi- 


nação linear de funções degrau unitário da seguinte maneira: 


e Para cada intervalo finito [a,b] onde f possui uma expressão diferente, multiplique a 


expressão de f pela diferença das funções degraus U(t — a) — U(t — b); 


e Para expressão definida num intervalo infinito [c, oo] somente multiplicar a expressão por 


Ult-c); 
e f(t) é a somatório de todos os produtos. 
Vamos ilustrar com os seguintes exemplos: 


Exemplo 2.7.11. Expresse 


Petro sensi<)> 
HD=4 é, se2<t<s5 


2tsent, set>5 


como uma combinação linear de funções de grau unitário. 


Solução: Observe que f(t) possui três expressões diferentes nos intervalos [0, 2), [2,5), e [5,00), 


logo temos duas diferenças, U(t) — U(t — 2) e U(t — 2) — U(t — 5) e três produtos, 
e(2-1+3) (2U(8) — Ult — 2)) para (f? — t+ 3) definida em [0,2); 
ec (Ut -D- Ult -— 5)) para é definida em [2,5); 


e 2tsent'U(t -5) para 2t sen t definida em [5, 00). 
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Portanto f(t) é a soma deste três produtos, 


ft) (P-++3) (UM) -U(t-2)) + (Ut — 2) - Ult— 5)) + 2º sen EU — 5) 


(P-++UMD+(e+P-t+BU(L-2)+(2tsent- e) U(L—5). 


Exemplo 2.7.12. Expresse 
Pe, se 0O<t<g3 
ft) =4 cost, se3<t<6 
2t, set>6 


como uma combinação linear de funções degrau unitário. 


Solução: Observe que f(t) possui três expressões diferentes nos intervalos [0,3), [3,6), e [6,00), 


logo temos duas diferenças, U(t) — U(t -3) e U(t —- 3) — U(t — 6) e portanto 


ft) Pe (U(L) — Ult — 3)) + cos t(U(t— 3) — U(t— 6)) + 2H U(t— 6) 


= PelUbD+(cost- Pe) U(t-3) + (2t- cos b) U(t- 6). 


2.7.5 Propriedade da Transformada de Laplace e Sua Inversa - Se- 


gundo Teorema da Translação 


Lema 2.7.6 ( Segundo Teorema da Translação). 


(a) Se Lf(t)) = F(s) ea for um número real qualquer, então, 


Liflt-)Ult-a) =eSF(s) ou LF UL-a)) = SLU(t+a)). 


(b) Se L!(F(s)) = f(t) então 





Demonstração. 
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OO 


e. f(t-a)Ult-a) dt 


Mt a) Ult — a) 


a 


e Ft a) Ult — a) dt+ E eft-a) Ult-a) dt 


et Ht— a) dt 


OO 


esto) F(u) du, (pela mudança u=t-a) 


SS 


Eee a e” f(u) du 


= pg CE(). 


(b) Seguir diretamente de (a). 


Exemplo 2.7.13. Calcule a Transformada de Laplace das seguintes funções 


cos (t —- 27/3), set>27/3 
(a) f(t) = sen tU(t - m) (b) f(t) = | 
0, set<21/3 
2, se0<t<2 
RP 
()HBD=4 —1, seZ2<t<3 (d) f(t) = | 
-1, t>5 
0, set>3 
Solução: 





(a) LlsentU(t-m)j=e“Llsen(t+m)|)=eTFLl-cost)=e”S. = = EE 
(b) Escrevendo f(t) em termos de funções degrau unitário temos que 
ft) = cos (t- 2/3) Ult - 2n/3). 
Portanto 
LI) = Llcos (t- 2n/3) Ut —2m/3)) = e PI cost) = e E 
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(c) Escrevendo f(t) em termos de funções degrau unitário temos que 
f(B = AU) — Ult 2)) = (Ult — 2) — Ult — 3) = 2U(8) - 3U( — 2) + UM 3). 
Portanto 


Lt) = LOUD) - LBU( — 2) + LU(t — 3)) 
1 es e 
pa pe 
S s s 








(d) Escrevendo f(t) em termos de funções degrau unitário temos que 
fe = (Ut) - Ult-5))-1U(t—5) = e! UM) — (6! +) U(t— 5). 
Portanto 


LD = Let UB) - Le! + 1) Ut 5) 


SS) = rio ae) 





= 1 —5s a 
- q (Lie e 2a] 
1 ese call 
o — — + — 
s+1 E s+1 | 
1 95 ES 
- se at CE 


Exemplo 2.7.14. Calcule 


See Ep (pe Sj —— a a | 
(es SS ua E) (en = re (s-52+16)' 








(a) £1 SS) =£Lfe*F(s)), onde F(s) = — > ft) = +. 


S— 
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Portanto 
q E = fe BE 
ea” te 2F(s)) 
= ft-DUL-2) 
= Dq). 
b) £! cia = Le Tip det == f) = 3t 
Lo | Le TPEO), onde F(s) = 5 = f() = cos (88) 
Portanto 
—Tis/2 
“(Ss = Lite ms2r(s)] 
= f(t-n/2) Ult- n/2) 
= cos3t— 7/2) Ult — /2) 
= — sen (3t) Ult — n/2). 
Ea seo =p Ads dalj=-s )= 3t 
(Lo [= Le ?sPF(o), onde F(9) = = f() = cos (3. 
Portanto 


e2rs/3 
E (5 + =) = Le 2sbr(s) 
= f(t- 2ms/3) Ult — 215/3) 
= cost — 275/3) U(t — 21715/3) 


= cos (3t) U(t — 275/3). 





= ER o pa o o 
(d) £ e EE =) =£L fe G(s)). onde G(s) = ESTE 
Seja F(s) = = então f(t) = L!F(s)) = cento 
Agora 
1 
G(s) = ESA =F(s-5). 
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Portanto, 
db) = L" no e a (F(s — 5)) = E f(t) = cam (48) 
(s—5)2 +42 4 
(pelo 1º teorema da translação). 
Portanto 


—3s 
= gt-3)U(t-3) 


1 
= qU(t —-3) 43) sen 4(t — 3). 


2.7.6 Derivada da Transformada de Laplace 


Lema 2.7.7. Se F(s) = Llf(t)) então 


Mais geral, 


Portanto, 





Exemplo 2.7.15. Determine a transformada de Laplace de É sen (5t). 


Solução: Temos que 
5 


Lf sen (5t)) = F(s) = DE 
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logo 


LIP sen (5t)) 


Il 


«Br 5 
qe [a E E] 
Ea —2s | 

ds |(s2 + 25) 
5(682 — 50) 
(522 +25)3 


II 





2.7.7 Resolução de Equações diferenciais pela Transformada de La- 


place 


2.7.7.1 Transformada de Laplace para Derivadas de uma função 


Lema 2.7.8. Seja y(t) uma função n- vezes derivável e satisfazem as hipóteses do Teorema 2.7.1, 


então 


d” 
L o E EO sp O) = 4H), 


onde Y(s) = Lly(t)). Portanto, 





Demonstração. Como F(s) := Lly(t)] = | e“! y(t) dt, temos que 
0 


LW(O) = / Esty(8) dt 


[es v(b| R +s / e“ty(t) dt integração por partes 
0 

0 — y(O) + sY(9) 

sY(s) — y(0) 


Il 
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Ly ()) = Lg(t)) onde g(t) = y'(t) 
= sG(s) — g(0) 
= sLlg(t — g(0) 
= sLy(b) — y'(0) 
= sisY(s) — y(0)) — y'(0) 
= sÊY(s) — sy(0) — y'(0). 


2.7.7.2 Resolução de equações diferenciais ordinárias de coeficientes constantes 


O problema de valor inicial da equação diferencial linear de ordem n : 


d"y det; dy 
E + EE ++ MT + any = g(t) 








com 
y(to) = yo 


dy ; 
“gr (to) =Y 


dy 
dtr—1 


pode ser resolvida através da transformada de Laplace como indicado no diagrama abaixo. 


(to) = no é 





Equação Diferencial ----- ==... 200020 > Solução da, . 
Equação Diferencial 


Transfbrmada Transformada 
de Laplace de Laplade Inversa 


Equação Algebrica O Soluçãoda 
Equação Algebrica 


Exemplo 2.7.16. Use a transformada de Laplace para encontrar y(t) : 


(a)y +4y=e E, W0)=2, 
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(b)y -y=1+té, y(0)=0, 

(Dy"+2y+y=0, v(O)=1, y(0)=1, 

(dy — 4y + 4y = el, y(0)=0, y'(0) =0. 
Solução: 


(a) y +4y=6*, y(O)=2. 


Tomando a transformada de Laplace: 


Lly)+4Lty) = Le) 
1 
=> sY(s) —- y(0) +4Y(s) = E 
> sY(s)-2+4Y(s) = — 
> (s+4)Y(s) = — +2 Equação Algébrica 


Resolvendo a equação algébrica por Y(s) temos que 








Tomando a Transformada de Inversa de Laplace, temos que 


cb ada 1 af do 
ví e le) t2L E 


> y(t) = fe po 





by -y=1+te, y(0)=0. 


Tomando a transformada de Laplace: 











> sY(s)-y(0)- Y(s) = : HE E — 
> sY(s)-Y(s) = : + E — 
> (s— DY(s) = : + E — 
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Resolvendo a equação algébrica por Y(s) temos que 














Tomando a Transformada de Inversa de Laplace, temos que 


ot e tes) 


1 

t 2 + 

e—-l+-—te. 
2 





«e 

E ear 
ja, 

e” 
| 


U 
a 
E 

HI 


()y'+2y +y=0, Y(0)=1, y(0)=1. 


Tomando a transformada de Laplace: 


Il 
O 


Ly) +2Lly) + Ly) 
= s2Y(5) = sy(O) — y (0) +2[sY(s)- y(O)) +Y(s) = O 


=> sY()-s-1+25Y(9)-2+Y(5) 


Il 
o 


= (8º +25 + 1D)Y(s) s+3 


Resolvendo a equação algébrica por Y(s) temos que 


o 8+3 E ss x 2. 
“242841 (s+12 s+1 (s+172 








Y(s) 


Tomando a Transformada de Inversa de Laplace, temos que 


A sf 
vi) = en) É (e) 
t 


= y(t) = ct +9te. 


(d) y” — 4y + 4y = Pet, y(O) = 0, y'(0) = 0. 
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Tomando a transformada de Laplace: 











Ly) -4Lly)+4Lty) = LP! 
= s*Y(s) — sy(0) — (0) — 4 [sY(s) — y(0)) + 4Y(s) = (s - 2)! 
g = Ec 
=> s?Y(s) — 4sY(s) + 4Y(s) = 6-2 
6 
=> (8º -45+4)Y(s) = 6-2 


Resolvendo a equação algébrica por Y(s) temos que 


6 E 
O DESA OCO (ER 


Tomando a Transformada de Inversa de Laplace, temos que 


sé —1 1 
“0 = is) 


6 
> y(t) = gre. 
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2.7.8 Tabela da Transformada de Laplace de Algumas Funções e 


suas propriedades 


= L!F(s)) 
ER) 
= a) 
no 


m 
= 
y 


I 
nd E 
gn+1 


a 





(a 


sen wt 


E + — 
É 
+w 
t cos wt 
t sen wt 
s 
Cao 
a 
pt 
O 
EE 
Lace ea 
— Ss 
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2.7.9 Exercícios 


1. Calcule a Transformada de Laplace £(f(t)], sendo: 





(a) f(t)=2e! —-5 sen (8t) + 7Ê. Resp: F(s) = — ico — + = 
(b) (1) = 2 cos (4t) + 4º *' sen (31). Resp: F(s) = ER E E 
(c) f(t) = 42e” + 6 cos (4t) — 10. Resp: F(s) = — - = - 
(d) f(t) = 5e! + 4 cos (6) — 2H. Resp: F(s) = — E E S 
360 2 9 


)=3Pe-4 )+9. Resp: F()=——— —- > + À 
(e) HB =3Pe sen (8t) +9 esp: F(s) 08 E é 


2. Escreva cada uma das funções em termos de funções degrau unitário. Ache a transfor- 


mada de Laplace da função dada. 














2» Dera 2 483 
(a) f(t) = Resp: f()=2-4U(t-3), F(s)=-—— 
=D bo s S 
1, O<t<4 
1 es e 5s 
(b) HD)=4 0, 4<t<5 Resp: HMD=1-U(t-9)+U(t-5), EB) = = E : 
1, t>5 
0, 0O<t<l 2 285 em 
E ; meio = = TSC 4 
(c) FC) = | É aa Resp: )=P-UL-1, )=D+ 5 
0, O<t<3n/2 = coBns/2 
(d) Ft) = s Resp: f(t) = sen t - U(t — Rs F(s) = —— 
sent, t>37/2 2 2+1 
t, O<t<2 te SRB 
(e) f(t) = Resp: (D=t-t-U(t-2), F(s)= 5" — 
0, t>2 s s s 
sent, O<t<27n 
(£) f(t) = Resp: f(t) = sen t-— sen t-Ul(t - 27), F(s) = 
0, t>27n 
1 Ens 


2+1 8241 
3. Calcule a inversa da Transformada de Laplace £ 1! (F(s)),, sendo: 


6 8 
decr 4. 
(eo Je) $$ s2+9 E s—5 
3 4s 5 

b) F()=-+>— 3 
Rojo 5 2+8l' 6+6P 


Resp: f(t) = 3P — É sen (3t) + 2e* 
E 3 





Resp: f(t) =3 +4 cos (9t) + Pes 
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S 


1 2 1 

(c) F(s) = E =aGEl Resp: f(t) = = cosh(2t) + 3 sinh(2t) + qe 
(d) F(9) = = E  Repufo=s eos 

4 2 6 2 
(e) F(s) = IDE: Resp: f(t) = -s Cos PES = Sen t+ Re 
(£) F(s) — e Resp: ft) — se sen (3t) 
(g) F(s) = a Resp: f(t) = ei cost 
(h) F(s) = — Resp: f(t) = é' cos (3) 





O) F)=—— Resp: fy)=e CI. uga-) 





Resp: HD=(1-t+D-U(t-2) 
4. Use a transformada de Laplace para resolver o problema de valor inicial dado. 


(a) y-3y=e!, y0O)=1 Resp: y(t) = 26! — e! 


4 
(b) y'-6y +9y=Pe!, yO)=2, y(0)=6 Resp: y(t) [+ 5)" 


() y' +16y= cos4t y0)=0,y(0)=1 Resp: y()= (+ uy den (id 
R 0, O<t<l 
(D) y +y=Ffb, vyO)=0, sendo f(t = 
5, t>1 
Resp: y(t) = [5 - 5e (DJ .U(t— 1) 
; 1, O<t<l 
(e) y +y=f(t, vy0)=0, sendo f(t) = 
-1, t>1 
Resp: ()=1-e!-2-U(t-1) +20 (CD. —1) 
) t, O<xt<l 
(D y +2y=f(t), y0)=0, sendo f(t) = 
0, t>1 
caio cha dradot Laguiono Lu ata de ED qt 
Resp: W()=-q+5!+7e q Ut-D-50-1)-UlL-D+ze U(t- 1) 


1, O<xt<l 
(9) y' +4y=f(tb, y0)=0, y(0)=-1, sendo f(t) = | 
0, t>1 


Resp: y(t) = : E : deste E E ; cos dli= 0) UE) : sen (28) 
(h) y” +4y= sent-U(t-27), yO)=1,y(0)=0 


Resp: y = cos (2t) — E sen 2(t — 270) - U(t — 27) + A sen (t — 27) - U(t — 277) 


MATB97 - Equações Diferenciais 


175 


0, O<t<rn 
Dy +y=ft, yvO)=0,y(0)=1, sendof(b)=41, n<t<2n 


0, t>2n 
Resp: y= sent +[1 — cos (t- 1)]-Ul(t - m) — [1 — cos (t — 27)] «Ult — 2m) 
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